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Abstract. We extend the perturbation theory of Visik, Ljusternik and Lidskii 
for eigenvalues of matrices, using methods of min-plus algebra. We show that 
the asymptotics of the eigenvalues of a perturbed matrix is governed by cer- 
tain discrete optimisation problems, from which we derive new perturbation 
formulas, extending the classical ones and solving cases which were singular 
in previous approaches. Our results include general weak majorisation in- 
equalities, relating leading exponents of eigenvalues of perturbed matrices and 
min-plus analogues of eigenvalues. 

1. Introduction 

Let Ae denote a, n x n matrix whose entries, which are continuous functions of 
a parameter e > 0, satisfy 

(1) (A)«j = a,je^'^' +o(e^'0 

when e goes to 0, where G C, and Aij G R U l+cxa}. (When Aij = +cxd, this 
means by convention that {Ae)ij is identically zero.) The goal of this paper is to 
give first order asymptotics 

CI ^ Ke^^ , 

with Ai G C \ {0} and A.^ G K, for each of the eigenvalues £g , . . . , of 

Computing the asymptotics of spectral elements is a central problem of pertur- 
bation theory, see |Kat95j and |Bau85j . For instance, when the entries of Ae have 
Taylor (or, more generally, Puiscux) series expansions in e, the eigenvalues have 
Puiseux series expansions in e, which can be computed by applying the Ncwton- 
Puiseux algorithm to the characteristic polynomial of Ae- The leading exponents 
Ai of the eigenvalues of A^ are the slopes of the associated Newton polygon: the 
difficulty is to determine these slopes from Ae- 
The case of a linear perturbation of degree one 

Ae = Ao + eb , Ao,bE C"^" , 

has been particularly studied. It suffices to consider the case where is nilpotent, 
which is the object of a theory initiated by Visik and Ljusternik |VL60j and com- 
pleted by Lidskii |Lid65j . Their result shows that for generic values of the entries 
of b, the exponents A; are the inverses of the dimensions of the Jordan blocks of 
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Ao- Then, the coefBcients Xi can be obtained from the eigenvalues of certain Schur 
complements built from the matrices and b. 

However, the construction of Visik, Ljustcrnik and Lidskii has many singular 
cases, in which the Schur complements do not exist, and so, their approach does 
not apply to non-generic situations, such as the case when the matrix b has a sparse 
or structured pattern. 

The problem of generalising the theorem of Visik, Ljustcrnik and Lidskii, i.e., of 
"categorising all possible behaviours as a function of the perturbation 5" , to quote 
the introduction of the article of Ma and Edelman |ME98| , has received much atten- 
tion. Their article solves cases where Ao and b have certain Jordan and Hessenberg 
structures, respectively. This problem is also considered in the survey of Moro, 
Burke, and Overton |MB097j . which includes a slight refinement of Lidskii's result 
together with an extension in special cases. Similar problems have been raised for 
matrix pencils, see in particular Najman |Naj99| . See also Edelman, Elmroth and 
Kagstrom |EEK97[ IEEK99j for a geometric point of view. Numerical motivation 
can be found there, as well as in the theory of pscudospectra, see Trefethen and 
Embree [TEOSj for an overview. In |Mur90j . Murota gave an algorithm to compute 
the Puiseux series expansions of the eigenvalues of a matrix whose entries are given 
by polynomials (or even formal series) in some indeterminate. This algorithm sheds 
some light on the problem raised by Ma and Edelman (although this problem is 
not considered there). 

In this paper, we use min-plus algebra to give elements of answer to the problem 
raised by Ma and Edelman. 

To describe our results, let us recall that the min-plus semiring, IR,nin, is the set 
MU {-|-cx)}, equipped with the addition (a, 6) i-^ min(a, 6) and the multiplication 
(a, b) a + b. Many of the classical algebraic constructions have interesting min- 
plus analogues. In particular, the characteristic polynomial function of a matrix 
B G ^n^n ' already introduced by Cuninghame-Green (CG83. , can be defined as the 
function which associates to a scalar x the permanent, in the min-plus sense, of the 
matrix xl ® B, where / is the min-plus identity matrix, "®" denotes the min-plus 
addition, and the concatenation denotes the min-plus multiplication. The perma- 
nent, in the min-plus sense, of a matrix B, is the value of an optimal assignment in 
the weighted bipartite graph associated with B. A result of Cuninghame-Green and 
Meijer CGM SQI shows that a min-plus polynomial function p{x) can be factored 
uniquely as p{x) — a{x © xi) ■ • ■ (x © Xn), where a,xi, . . . ,Xn S Rmin- The num- 
bers Xl, . . . , Xn, which coincide with the points of non-differentiability of p (counted 
with appropriate multiplicities), are called the roots or corners of p. The sequence 
of roots of the min-plus characteristic polynomial of a matrix B G MJJ,^," can be 
computed in polynomial time, by solving 0(n) optimal assignment problems, as 
shown by Burkard and Butkovic [BB03l| . The reader seeking information on the 
min-plus sem i ring may consult pT79l IMS92I |BCUQ92| IMax94l l(X;95l I(;un98l 
IK M97I lGP97l IPin98l KWMl \^M0^ . 

We assume that Ae is given by This allows one to handle the case of a 
perturbed matrix Ae ~ Aq + eb, where the matrix b is non-generic. 

The first main result of the present paper. Theorem 13. 81 shows that the sequence 
of leading exponents of the eigenvalues of the matrix Ae is weakly (super) majorised 
by the sequence of roots of the min-plus characteristic polynomial of the matrix 
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of leading exponents of Ae, and that the equahty holds for generic values of the 
coefficients Uij. 

The proof of Theorem 13.81 relies on a variant of the Newton-Puiseux theo- 
rem in which the data arc only assumed to have first order asymptotics, that 
we state as Theorem 13.11 in a way which illuminates the role of niin-plus alge- 
bra. We consider the branches y{e) solutions of the equation V{e,y{e)) — 0, 
where 7^(e, 1^) = X]J=o (^)^"' ^^"^ ^^'^ ^j (^) ^^'"^ continuous functions, such that 
Vjie) = pje^i + o(e^j ), with pj G C and G K U {+cx}}. Wc characterise the 
cases where this information is enough to determine the first order asymptotics of 
the branches 3^i(e), . . . , yn(e). Then, the leading exponents of the branches are 
precisely the roots of the min-plus polynomial P{Y) = ^"^QPjY-'- the leading 
exponents of the classical roots are the min-plus roots. Note that in this case, by 
Legendre-Fenchel duality, the roots of the min-plus polynomial P{Y) are precisely 
the slopes of the Newton-Polygon classically associated to V{e,Y). Note also that 
the generic equality in Theorem 13.81 could be derived from Murota's combinatorial 
relaxation technique |Mur90j . which uses a parametric assignment problem, whose 
value is exactly the min-plus characteristic polynomial. 

Theorem 13.81 determines the generic leading exponents of the eigenvalues of Ae, 
but it does not determine the coefficients A^. To compute these coefiicients, we 
define, in terms of eigenvalues of min-plus Schur complements, a sequence of critical 
values of A, that we characterise as generalised circuit means. We show that the 
sequence of roots of the min-plus characteristic polynomial of A is weakly majorised 
by the sequence of critical values of A (Theorem I4.6|l , and we characterise the 
equality case in terms of the existence of disjoint circuit covers, or perfect matchings, 
in certain graphs. 

Our second main result, Theorem 15. II shows that, in the equality case of Theo- 
rem l4.6l the coefficients A; can be obtained in terms of eigenvalues of certain Schur 
complements constructed from the matrix a. The theorem of Visik, Ljustcrnik and 
Lidskii is a special case of this result fCorollarv l7.1|l . Wc give in Section [731 ex- 
amples of singular cases which can be solved by Theorem 15.11 In the remaining 
singular cases, different methods should be used, along the lines of jMurQOI rABG04j . 

We also prove an asymptotic result for eigenvectors. Theorem 16. II which is anal- 
ogous to Theorem 15.11 However, the combinatorial characterisation of the cases 
where Theorem 16.11 determines the generic asymptotics of all the entries of eigen- 
vectors is lacking, see Section Note that even when the first order asymptotics 
of an eigenvalue is determined, a detailed asymptotic information on Ae may be 
needed to determine the first order asymptotics of the corresponding eigenvector, 
as shown in our earlier work |ABG98j which concerns the special case of the Perron 
eigenvector. 

The present results provide a new illustration of the role of min-plus algebra in as- 
ymptotic analysis, which was recognised by Maslov |Mas73l Ch. VIII]. He observed 
that WKB-type or large deviation type asymptotics lead to limiting cc|uations, like 
Hamilton- Jacobi equations, satisfying some idempotent superposition principle. So, 
min-plus algebra arises as the limit of a deformation of usual algebra. This obser- 
vation is at the origin of idempotent analysis ;MS92', 'Klvig? , 'LMS01\ It has been 
used by Dobrokhotov, Kolokoltsov, and Maslov | .DKM92. ..KM97j to obtain precise 
large deviation asymptotics concerning the Green kernel and the first eigenvalues of 
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a class of linear partial differential equations, with application to the Schrodinger 
equation. 

The same deformation has been identified by Viro jVirOl) . in relation with the 
patchworking method he developed for real algebraic curves. It appears in several 
recent works in "tropical geometry", in particular, by Mikhalkin [MikOll fMik03j . 
Forsberg, Passare, and Tsikh [FPTOOj . Passare and RuUgard )PR,04j . and Speyer 
and Sturmfels |SS04j , following the introduction of amoebas of algebraic varieties by 
Gelfand, Kapranov, and Zelevinsky |GlKZ94j . In these works, the relation between 
Newton polytopes and min-plus or "tropical" polynomials is apparent. We use the 
same relation in the version of the Newton-Puiseux theorem concerning first order 
asymptotics that we stated as Theorem II 

Relations between max-plus algebra and asymptotic problems have also appeared 
in other contexts. Puhalski jPuhOlj applied idempotcnt techniques to large de- 
viations theory. Friedland |Fri86j observed that the max-plus eigenvalue can be 
obtained as a limit of the Perron root. Olsdcr and Roos |OR88j and De Schut- 
ter and Dc Moor |DSDM98) used asymptotics theorems to derive certain max-plus 
algebraic identities. 

Finally, wc note that Theorem 15 . 1 1 was announced in |ABG!01j . and that the role 
of min-plus roots in the Newton-Puiseux theorem was mentioned in |GP01| . 

2. Preliminaries 

In this section, we recall some classical facts of min-plus algebra and show pre- 
liminary results. See for instance |BCOQ92| for more details. 

The min-plus semiring, Mmin, is the set R U {-|-cx)} equipped with the addition 
{a,b) ^ a (B b ~ min(a, 6) and the multiplication (a, 6) >-^a^b^a + b. We shall 
denote by (D = -l-oo and 1 = the zero and unit elements of Mmin, respectively. We 
shall use the familiar algebraic conventions, in the min-plus context. For instance, 
if A, B are matrices of compatible dimensions with entries in Kmin, {AB)ij = [A® 
B)i3 = 0fc AikBkj = mmk{A,k+Bkj), A^ = etc. Moreover, if x S Rmi„\{(D}, 

then is the inverse of x for the ® law, that is — x, with the conventional 
notation. We shall also denote by Mmin the complete min-plus semiring, which is 
the set RU{±oo} equipped, as Mmin, with the min and -I- laws, with the convention 
+00 + {—oo) = — oo -I- (+oo) = -l-oo. 

2.1. Min-plus spectral theorem. To any nxn matrix A with entries in a semir- 
ing S, we associate the directed graph G{A), which has nodes 1, . . . , n and an arc 
(i, j) \i Aij ^ (D, where (D denotes the zero element of S. We say that A is irreducible 
if G{A) is strongly connected. 

We next recall some results of min-plus spectral theory: the min-plus version 
of the Perron-Frobenius theorem has been discovered, rediscovered, precised or 
extended, by many authors |(X;79llXAl767llRom67linM77l|CDQV88[IMS9^^ Recent 
presentations can be found in |IjCUg92| IHnoHl l(;P97l |ljap98| lAciw] . 

Theorem 2.1 (Min-plus eigenvalue, see e.g. |BCOQ92| Th. 3.23]). An irreducible 
matrix A G (Kmin)"'^" has a unique eigenvalue: 



(2) 



n 
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With the usual notations, (0 can be rewritten as 

^n»2 + ' ■ ' + Ai 



-^{A) = min min 



l<k<n ii,...,ik k 

lip = (io,«i, . . . ,«fe) is a path ofG{A), we denote by \p\a = ^loii ^ ^^ik-iik the 

weight of p, and by |p| = fc its length. Since any circuit of G{A) can be decomposed 
in elementary circuits, which are of length at most n, /Omin(^) is the minimal circuit 
mean: 

(3) Pmin{A) = min . 

c circuit in G(A) |C| 

We say that a circuit c = (ii, «2, • ■ • , ^fc, ^i) of G{A) is critical if c attains the 
minimum in Q, and we call critical the nodes and arcs of this circuit. The critical 
nodes and critical arcs form the critical graph, G'^(A). We call critical classes the 
strongly connected components of ^{A). Wc will also use the name "critical class" 
for the set of nodes of a critical class. 

The Kleene 's star of a matrix A G K^^^jl* is defined by 

A*=I®A®A'^®---e 1",^„" , 

i.e. {A*)ij — inffe>o(^'^)y , where / = A'' is the identity matrix (we shall use the 
same notation / for the identity matrix of K^^j^*, and for the identity matrix of 
i^nxn^ for any n). 

Proposition 2.2 (See e.g. |BCOQ92| Th. 3.20]). All the entries of A* are > -oo 

if, and only if, PminiA) > 0. Moreover, when PminiA) > 0, 

A* = / ® A © • • • © A"~^ . 

Theorem 2.3 (Min-plus eigenvectors, see e.g^ |BCUQ92| Th. 3.100]). Let A e 

^min '^'"^ irreducible matrix, and consider A = pnii,i(A)~^A. Any eigenvector 
of A is a linear combination of the columns A* ^ corresponding to critical nodes j. 
More precisely, if we select (arbitrarily) one node j per critical class and take the 
corresponding column A* p we obtain a minimal generating set of the eigenspace of 
A. 

(In Theorem 12.31 and in the sequel, we write A* ^ the j-th column of {A)* .) 

Given a matrix A £ l^mhT ^^"^ ^ vector V S ^JJ^jj^, we define the saturation 
graph, Sat(A, V), which has nodes \, . . . ,n, and an arc («, j) if {AV)i = AijVj (that 
is {AV)i = Aij + Vj with the usual notations). The following simple result relates 
the critical graph and the saturation graph. 

Proposition 2.4 (See e.g. |BCUQ92| Th. 3.98]). Let A e M;\^„" be an irreducible 
matrix with eigenvalue a, and let V S R^in \ {"^l- V ~ o;^,' then the strongly 
connected components of Sat(A, V) are exactly the strongly connected components 
ofG^A). 

In fact, Theorem 3.98 of |BCOQ92| only shows that any circuit of the saturation 
graph belongs to the critical graph, but the converse is straightforward. 

The following elementary result is a special version of a maximum princi- 
ple for ergodic control problems, see |AG!03I Lemma 3.3] for more background, 
and jCTGG!99l Lemma 1.4] for a proof in the min-plus case. 
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Proposition 2.5. Let A G l^mii" irreducible matrix with eigenvalue a, and 

let V e If AV > aV, then {AV)i = aVi for all critical nodes i of A. 

The saturation graphs associated to the generators of the eigcnspace have a 
remarkable structure. Say that a strongly connected component C of a graph is 
final if for each node i, there is a path from i to C, and if there is no arc leaving C. 

Proposition 2.6. Let A e Kmh" '^'^ irreducible matrix with eigenvalue a, let 
A = a^^A, let C be a critical class of A, and let V be an eigenvector of A. The 
following assertions are equivalent: 

(1) V is proportional to A*p for some j G C; 

(2) C is the unique final class ofSat{A,V). 

Proof. We first prove d =^ |21 It is enough to consider the case when V = A* . 
Since A is irreducible, all the entries of A* are < +oo. Moreover, since jOmin(^) = 0, 
Proposition E21 yields A* = J © A © ■ • ■ © A""^ Hence, for aU i ^ j, there exists 
a path p — (io — i,ii, . . . ,ik — j) from i to j, with length 1 < fc < n — 1, and 
minimal weight, that is A*j = A^yi-^ ■ ■ ■ Ai^^ -^i^. By Bellman's optimality principle, 
for all < / < TO < fc, the sub-path {ii, . . . ,im) has minimal weight: A*^^^ ~ 
• ■ • Then, 1*^- = 1,;,,;,^^ A* ^^j-, that is, aV^, = A,,^,,^^^V^,^,, and 

{ii, ii+i) £ Sat (A, V) for alH = 0, . . . , — 1. So for i ^ j, there is a path from i to 
j in Sat{A,V). 

Assume, by contradiction, that there exists k ^ C and I ^ C such that (A:, I) £ 
Sat(A, V). Since I ^ j, there is a path from I to j in Sat(yl, V), and since C is a 
strongly connected component of Sat{A,V) (by Proposition 12.4(1 . there is a path 
from j to k in Sat(A, y), which yields a circuit of Sa,t{A,V) passing through C 
and k ^ C . This contradicts the fact that C is a strongly connected component of 
Sat(A, V). 

We finally prove[5]==>n Assume that C is the unique final class of Sat (A, V), 
and let us fix j G C . Then, for each i, we can find a path (iq = i, . . . ,ik = j) 
from i to j in Sat(A, y), so that Vt„ = A*^^^Vi^, V^^_-, = A*^_^^^^V^^. Hence, 
Vi = A*^^^ ■ ■ ■ A*^_^i^Vj < A*jVj. The other inequality holds, since V = AV implies 
V = A*V. Thus, V = 1* V, is proportional to 1* . □ 

2.2. Min-plus polynomials. We recall here some results about formal polyno- 
mials and polynomial functions over Mmin, and in particular a min-plus analogue 
of "the fundamental theorem of algebra" , which is due to Cuninghame- Green and 
Meijer ICXtMSO] . The connection between the min-plus evaluation morphism and 
the Fenchel transform, was already observed in |CGNQ89| and |BCOQ92| Sec- 
tion 3.3.1]. 

We denote by Mniin[Y] the semiring of formal polynomials with coefficients in 
Rmin hi the indeterminate Y: a formal polynomial P G R,„i,i[Y] is nothing but a 
sequence {Pk)k&i S l^min such that = (D for all but finitely many values of k. 
Formal polynomials are equipped with the entry-wise sum, (P © Q)^ = Pk ® Qk, 
and the Cauchy product, {PQ)k ~ ®o<i<fe PiQk-i- As usual, we denote a formal 
polynomial P as a formal sum, P = ©/^o-^fc^*^- define the degree and 

valuation of P: degP = sup{fc G N | P^ 7^ C}, valP = inf{fc G N | Pfe 7^ C} 
(degP — —00 and valP = +(X) if P = 0). To any P G Mmin[Y], we associate the 
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polynomial function P : Mmin ^ Rmin, y ^ P{y) ~ ®^=o -P^y'^ , that is, with the 
usual notation: 

(4) P{y) = min(Pfc + ky) . 

We denote by Mmin{Y} the semiring of polynomial functions P. Contrary to 
the case of real or complex polynomials, the evaluation morphism, Mniin[Y] — > 
Kmin{Y}, P 1-^ P is not injcctivc. Indeed, the evaluation morphism is essentially a 
specialisation of the Fenchel transform over R: 

^ : r"* ^ r", Hf)iy) = ^M^y - fi^)) , 



since, for all y G R, P{y) = — JF(P)(— y), where P is extended to a function 
(5) P:R^R, a;^P(x), withP(a;) = 



Pfe ifx^keN , 
+00 otherwise 



It follows from Q that P is a concave nondecreasing function with integer slopes. 

In the sequel, we denote by vex/ the convex hull of a map / : R R, and we 
denote by P the formal polynomial whose sequence of coefficients is obtained by 
restricting to N the convex hull of the map P : R ^ R. Thus, Pk = (vexP)(fc). The 
following result is a special case of the Legendre-Fenchel inversion theorem |Roc70l 
Section 12]. 

Proposition 2.7. If P £ R,nin[Y], then P is the minimal formal polynomial Q 
such that Q ~ P, we have P = P, and P is given by 

Pk=sup{-ky + P{y)) . 

Theorem 2.8 ( |BCOQ92| Th. 3.43, 1 and 2]). A formal polynomial of degree n, 
P G Rmin [Y] , satisfies P = P if and only if there exist Ci < ■ ■ ■ < Cn G Rmin such 
that 

P = Pn{y®Ci)---{Y(BCn) . 

The Ci are unique and given, by: 

,c\ \ Pn-i{Pn-i+l) ^ if -fri-i+1 7^ (D , 

(6) = < . for I = 1, . . . ,n . 

I (U otherwise, 

The min-plus analogue of the fundamental theorem of algebra due to 
Cuninghamc- Green and Meijer can be obtained by applying Theorem 12.81 to P, 

since P ~ P and P = P. 

Theorem 2.9 f |CGM80j ). Any polynomial function P e Rmin{Y} can be factored 
in a unique way as 

(7) P{y) ^ Pn{y ® ci) ■ ■ ■ {y ® Cr,) , 
with ci < • • • < c„ . 



The Ci are called the roots of P. (In |(XtM80) . the term corners is used 
as a synonym of root, we use the term of root which makes the analogy with 
classical algebra clearer.) The multiplicity of the root c is the cardinality of the 
set {j e {l,...,n} I Cj — c}. We shall denote by R(P) the sequence of roots: 
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R(P) = (ci, . . . , Cn). By extension, if P G Minin[Y] is a formal polynomial, we will 
call roots of P the roots of P: R(P) = R(P). By Proposition O R(P) = R(P). 
Geometrically, the function P is the restriction to N of the convex function vexP, 
which is piccewisc affinc on its support, [valP, degP], and P is concave, piecewise 
affine. 

Proposition 2.10. The roots c (IzR of a formal polynomial P G IR„iin[Y] are exactly 
the points at which P is not differentiable. They coincide with the opposites of the 
slopes of the affine parts ofveyiP : [valP, degP] — > R. The multiplicity of a root 
c G R is equal to the variation of slope of P at c, P'{c^) — P'{c~^), and it coincides 
with the length of the interval where vcxP has slope — c. Moreover, (D is a root 
of P if and only if P'((D~) := limc^+oo -P'(c) ^ 0. In that case P'((D") is the 
multiplicity of (D, and it coincides with val P. 

Proof. The characterisation of the roots and of their multiplicities in terms of P 
is due to Cuninghame- Green and Meijer | CGM80| . It can be deduced from lO, 
since when c G R, Qiy) {y © c)*"' = fcmin(j/,c) has c as unique point of non 
differentiability, with Q'{c^) = k and Q'{c'^) = 0. The case where c = (D is 
a straightforward consequence of ((T)). The characterisation of the roots and of 
their multiplicities in terms of vexP follows from 0, since when Ci G R, = 
P„_i — P„_i+i = (vcxP)'(a;) for all x G (n — i, n — i + 1); and Ci = (D =^ Pn~i = 
P„ci • • • Ci = (D. □ 

The duality between roots and slopes in Proposition 12.101 is a special case of 
the Legendre-Fenchel duality formula for subdifferentials: — c G d{vexP){x) <J=> 
X G dJ-{P){—c) X & d~^P{c) where d and d'^ denote the subdifferential and 
superdifferential, respectively |Roc70l Th. 23.5]. 

Lemma 2.11. Let P ~ ©"=0 ^i^^ € Rinin[Y] he a formal polynomial of degree n. 
Then, R(P) = (ci < • • • < c„) if, and only if P > P„(Y ci) • • • (Y c„) and 

(8) Pn-i = P„ci • • • Ci for all i G {0, ?i} U {i G {1, . . . , rt - 1} | a < Ci+i} . 

In particular, Pn-i = Pn-i holds for all i as in ((HJ. 

Proof Wc first prove the "only if" part. If R(P) = (ci < • • • < c„), then P = 
Pn(y ® ci) ■ • ■ (Y © c„) and Pn-i = PnCi ■ ■ ■ Ci for all i = 1, . . . n. Recall that 
P defines a map x ^ P{x) by ©. By definition of vcxP, the epigraph of vexP, 
epi vex P, is the convex hull of the epigraph of P, cpi P. By a classical result jR.ocTOI 
Cor 18.3.1], if S' is a set with convex hull C, any extreme point of C belongs to S. 
Let us apply this to 5' = epiP and C = epi vex P. Since Pn-i ~ PnCi ■ ■ - Ci, the 
piecewise affine map vexP changes its slope at any point n — i such that Ci < Ci+i. 
Thus, any point (n — i, vexP(n — i)) with Ci < Ci+i is an extreme point of epi vex P, 
which implies that (n — i,vexP(n — i)) G epiP, i.e., P„_i < vexP(n — i) = Pn-i- 
Since the other inequality is trivial by definition of the convex hull, wc have Pn-i = 
Pn-i- Obviously, P and P have the same degree, which is equal to n, and they 
have the same valuation, k. Then, (n, vexP(n)) and (fc,vexP(fc)) are extreme 
points of epi vex P, and by the preceding argument, P„ = P„, and Pk = Pk- 
Hence, Pq = Po, if A: = 0, and Pq = Po = +oo, if fc > 0. We have shown JHJ, 
together with the last statement of the lemma. Since P„ ~ P„ and P > P, wc also 
obtain P > P„(Y © ci) • • • (Y © c„). 
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For the "if" part, assume that P > P„(Y e ci) • • • (Y ® c„) and that (jHJ holds. 
Since Q = P„(Y Ci) • • • (Y ® c„) is convex, and the convex huU map P ^ P is 
monotone, we must have P > Q = Q. Hence, P > P > Q and since Pn-i = Qn-i 
for aU i as in ©, we must have Pn-i = Qn-i, thus vexP(n — i) = vex.Q(n — i) 
at these i. Since vexP is convex, since vexQ is piecewise afhne and vex(5(j) = 
vexP(j') for J at the boundary of the domain of vexQ and at all the j where 
vexQ changes of slope, we must have vexP = vexQ. Hence P = Q = Q and 
R(P)-R(P) = R(Q) = (ci,...,c„). □ 

The above notions are illustrated in Figure^ where we consider the formal min- 
plus polynomial P = Y"^ © 5Y^ © 6Y © 13. The map j t-^ Pj , together with the map 
vexP, are depicted at the left of the figure, whereas the polynomial function P is 
depicted at the right of the figure. We have P = Y^eSY^eGYelS = (Y©3)2(Y®7). 
Thus, the roots of P are 3 and 7, with respective multiplicities 2 and 1. The roots 
are visualised at the right of the figure, or alternatively, as the opposite of the slopes 
of the two line segments at the left of the figure. The multiplicities can be read 
either on the map P at the right of the figure (the variation of slope of P at points 
3 and 7 is 2 and 1, respectively), or on the map vexP at the left of the figure (as 
the respective horizontal widths of the two segments). 




Figure 1. The formal min-plus polynomial P = Y3©5Y2e6Y©13 
and its associated polynomial function P. 



2.3. Schur complements. We recall here the definitions of conventional and min- 
plus Schur complements. We shall consider matrices indexed by "abstract indices" : 
if L and M are finite sets and S* is a semiring, & Lx M matrix with values in S is an 
element A of S^^^^ and the entries of A are denoted by Aij with i € L and j G M. 
Moreover, for all J C P and K C M, we denote by Ajk the J x K submatrix 
of A: Ajk = {Ajk)ji£j^kGK- This definition applies to n x n matrices by taking 
L = M = {1, . . . ,n}. Graphs oi L x L matrices A are defined as for n x n matrices 
(see Section im with the only difference that the set of nodes is L. 

Definition 2.12. Let C C P be finite sets, and let = P \ C. If a is a P x L 

matrix with entries in C, and if acc is invertible, the Schur complement of C in a 
is defined by 

Schur(C, a) = unn - aNciaccY^ acN ■ 
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Definition 2.13. Let C C L be finite sets, and let N ^ L\C. If A is a L x L 
matrix with entries in Mmin, A G Mmin \ {©}, and Pmini^^^Acc) > 0, the min-plus 
X-Schur complement of C in A is defined by 

(9) Schur(C, A, A) = Ann © ANciX-'Acc)* X'^Acn ■ 

When A = 1 = 0, we shall simply write Schur(C, A) instead of Schur(C, 1, A). 

In fact, in the sequel, we shall mostly use min-plus Schur complement corre- 
sponding to A = Pmin{A). The goal of the insertion of the normalising factors in 
is to get the following homogeneity property: 

(10) Schur(C, fiX, pA) = fi Schur(C, A, A) , 

for all A, G M such that A < Pmin(^cc) a-nd pX < Pmin{Acc')- 

Using the same symbol, "Schur" , both for conventional and min-plus Schur com- 
plements is not ambiguous: considering min-plus Schur complements of complex 
matrices, or conventional Schur complements of min-plus matrices, would be mean- 
ingless. 

Both min-plus and conventional Schur complements satisfy 

(11) Schur(C U C", a) = Schur(C, Schur(C", a)) 

for all L X L matrices a, and for all disjoint subsets of indices C, C C L, provided 
that the Schur complements are well defined (if Schur(C",a) is well defined, then 
the left hand side of (|11|) exists if, and only if, its right hand side exists). Of 
course, (lllll is a classical Gaussian elimination identity, which is well known, both 
in conventional algebra and in the min-plus algebra (the left hand side and the 
right hand side of Hll|) are unambiguous rational expressions, with elementary 
interpretations in terms of paths, see for instance |Lal79j for more background). 

Finally, if K C L and if b is the K x K submatrix of a, we shall sometimes write 
abusively Schur(6,a), instead of Schur(i^, a). 

We now give some graph interpretations of the weights and eigenvalues of min- 
plus Schur complements. Let G be a graph with set of nodes L, let C be a subset 
of L and set N ^ L\C. For all paths p = (ig, . . . , ik) of G, we denote by \p\c the 
number of arcs of p with initial node in C, i.e., \p\c = #{0 < m < fc — 1 | im £ C}, 
where # denotes the cardinality of a set. (All the path interpretations below have 
dual versions, obtained by replacing "initial" by "final".) We also denote by pfl C 
the subsequence of p obtained by deleting the nodes not in C {pH C need not be 
a path of G). The following classical interpretation of Schur complements is an 
immediate consequence of the graph interpretation of the star. 

Lemma 2.14. Let C <Z L be finite sets, and let N ~ L\G . Let A be a L x L matrix 
with entries in Mmin, and X £ Mmin \ {©} be such that pminiAcc) ^ X. Then, p is 
a path in G(Schur(C, A, A)) if, and only if, there exists a path p' in G{A) with the 
same extremal nodes as p and such that p' H N = p. Moreover, for all paths p in 
G(Schur(C, A, A)), we have 

|p|schur(c,A,A) = min|p'|A - X\p'\c , 

where the minimum is taken over all the paths p' of G{A) that have the same 
extremal nodes as p and satisfy p' H N = p. In particular, c is a circuit in 
G(Schur(G, A, A)) if, and only if, there exists a circuit c' in G{A) such that 
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c' n N = c; and for all circuits c in G(Schur(C, A, A)), we have 

|c|schur(aA,A) = min \c'\a - X\c'\c , 

where the minimum is taken over all the circuits c' of G{A) such that c' Cl N ~ c. 

Proposition 2.15. Let C C L be finite sets, and let N = L\C. Let A be a L x L 
matrix with entries in Mmin, and X £ M,nin\{(D} be such that PminiAcc) > A. Then, 

(12) p„,in(Schur(C, A, A)) - min ^''[^f'^^^'^^ 

\c'\ ~ \c'\c 

where the minimum is taken over all the circuits c! of Gi^A) which are not included 
in C. Moreover, c is a critical circuit of Schm:{C, X, A) if, and only if, there exists 
a circuit d of G{A) such that c' D N = c and c' minimises (|12|l . 

Proof. Using (0) and Lemma [2. 141 we get 

Pmin(Schur(C, A, A)) = min I lschur(c,A,A) 

c circuit in C 

/ . \c'\a-X\cJ\c 
— mm mm 

c circuit in N \c' circuit of G{A),c'nN=c \c\ 

\c'\a-X\c'\c 

c' circuit oi G{A),c'nN^0 \c'\ — \c\c' 

since |c' n -/V| = |c'| — jc'lc for all circuits c'. This yields H12|) . If c is a criti- 
cal circuit of Schur(C, A, then /9niin(Schur(C, A, A)) = (|c|schur(c.A,A))/|c| and 
by Lemma 12.141 there exists a circuit c' of G{A) such that c' n N ~ c and 
c|schur(c,A,A) = W\a — X\c'\c- Siuce in that case, |c| = |c' n = |c'| — \c'\c, 
we deduce that c' minimises H12f) . Conversely, if c' minimises (I12II , then, c — c' CiN 
is nonempty and by Lemma [2. 141 c is a circuit of G'(Schur(C, A, A)). Moreover, by 
Lemma 12.141 again, 

/C U ir^ \ AW ^ l'^lschur(C,A,^) , |c'U-A|c'|c , . 

Pmin(Schur(C, A, A)) < p— '- < —— = pnun{Schur[C, X, A)) , 

\c\ \c'\ - \c'\c 

thus c is a critical circuit of Schur(C, A, A). □ 

Note that if c' is a circuit in C, that is if the denominator in H12() is zero, the 
numerator is necessarily nonnegative, since A < Pmin(^cc)- 



3. MiN-PLUS POLYNOMIALS, NeWTON-PUISEUX THEOREM AND GENERIC 
EXPONENTS OF EIGENVALUES 

3.1. Preliminaries on exponents and general assumptions. Let C denote the 
set of continuous functions / from some interval (0, eo) to C with eq > 0, such that 
< e"*^ on (0, eo), for some positive constant k. Since all the properties that we 
will prove in the sequel will hold on some neighbourhoods of 0, we shall rather use 
the ring of germs at of elements of C, which is obtained by quotienting C by the 
equivalence relation that identifies functions which coincide on a neighbourhood of 
0. This ring of germs will be also denoted by C. For any germ f G C, we shall 
abusively denote by /(e) or /<. the value at e of any representative of the germ /. 
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We shall make a similar abuse for vectors, matrices, polynomials whose coefRcients 
are germs. We call exponent of / G C: 

(13) e(/) liminf l^^LM g R y {+00} . 

£^0 loge 

Wc have, for all f,g^C and A G C, 

(14) e(/ + .g) >min(e(/),e(.9)) , 

(15) e(/.g) >e(/) + e{g) , 

with equality in (|14() if e(/) ^ e{g) and equality in H15|l if the liminf in the definition 
of e(/) or e{g) is a limit. Thus, / 1— > e(/) is "almost" a morphism C Rmin- In 
the sequel exponents will be considered as elements of Mmin, so that H15|) will be 
written as e{fg) > e{f)e{g). An element / G C is invertible if, and only if, e(/) ^ (D 
(or equivalently, if there exists a positive constant such that |/(e)| > e*^)- If / is 
invertible, its inverse is the map : e 1— > /(e)~^ and we have e(/^^) < e(/)~^ 
with equality if, and only if, the liminf in the definition of e(/) is a limit. 
We shall say that / G C has a first order asymptotics if 

(16) /(e) - ae^, when e ^ 0+ , 

with either A G M and a G C \ {0}, oi A = +00 and a G C. In the first case, (|1()|) 
means that lime_,o e~"^/(e) ~ a, in the second case, p6|l means that / = (in a 
neighbourhood of 0). We have: 

(17) /(e) ~ ae^ =^ e(/) = A , 

and the liminf in H13() is a limit. We shall also need an equivalence notion slightly 
weaker than ~. If / G C, a G C and A G Mmin, we write 

(18) /(e) ^ ae^ 

if /(e) = ae^ + o(e'^). If A G M, this means that lim^^o e~^/(e) = If ^ = +00, 
this means by convention that / = 0. If a or A = +00, then /(e) 2± ae^ if, and 
only if, /(e) ~ ae^ and in that case e(/) = A. In general, 

(19) /(e) ^ ae^ =^ e(/) > A . 

Conversely, e(/) > A /(e) ~ Oe^^. Of course, in (|18() . ae^ must be viewed 

as a formal expression, for the equivalence to be meaningful when a = and A G M. 
In (|17|) . however, ae^ can be viewed either as a formal expression or as an element 
of C. 

Throughout the paper, we consider a matrix A G C"^" and we shall assume that 
the entries {Ae)ij of Ae have asymptotics of the form: 

(20) {Ae)ij ^ ajje'^*^ for some matrix a = {aij) G C"^", 

and for some irreducible matrix A — {Aij) G K-mhT- 

(The case where A is reducible is a straightforward extension.) Under rather general 
circumstances (see Section IS^ . the eigenvalues C\,. . . oi A^ belong to C and 
have first order asymptotics: 

(21) C.^Ke^^ . 

We next relate the sequence (Ai,...,A„) with two sequences constructed by 
using only the information on the exponents of the entries {At)ij of the matrix Ae 
given by the Aij . 
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3.2. First order Newton-Puiseux theorem and min-plus polynomials. The 

usual way to compute the Ai in (|21(l is to use the classical Newton-Puiseux theorem. 
We state here a general first order version of this theorem in a way which illuminates 
the role of min-plus algebra. 

For any formal polynomial with coefficients in C, Vie, Y) = Sj=o (^)Y"' ^ , 
we define the min-plus polynomial of exponents: 

n 

The transformation of ordinary polynomials to min-plus (or "tropical" ) polynomial 
by the map e is instrumental in works on amoebas (for instance, a very similar 
definition is given in [SH^). 

Recall that to P = 6(7-") is associated the polynomial function P and the convex 
formal polynomial P, as in Section For instance, to = Y'^ + e^Y^ — e^Y -I- e^^ 
corresponds the formal min-plus polynomial P — e{P) = Y"^ © 5Y^ © 6Y -f 13 
represented in Figure ^ 

Theorem 3.1 (First order Newton-Puiseux theorem). Let V = j=o ■'^i (^)^"' ^ 
C[Y] such that Vn = 1- The following assertions are equivalent: 

(1) There exist 3^i,...,3^„ S C such that 3^i(e), . . . , y„(e) are the roots of 
Vie,!/) — counted with multiplicities, and 3^1,..., 3^n have first order 
asymptotics, 3^j(e) ~ Uj^^^ with Yi < ■ ■ ■ <Yn; 

(2) There exist p = J2"=oPj^^ ^ and P = ©^^o ^ IRmi„[Y] sat- 
isfying Vj{e) ~ Pji^' , j = 0,...,n, with p„ = 1, P„ = 1, po ^ or 
Pq = 0, and Pn-i ^ for all i G {1, . . . , n — 1} such that Ci < ci+i, where 
(ci < ••• < c„) = R(P). 

When these assertions hold, we have e(P) > P, e{V) = P. and R(e(P)) = R(P) = 
(ci < ••• < c„) = (Yi < ••• < Y„). Moreover, if c E M is a root of P with 
multiplicity k and Ci ~ ■■■ = Ci+k-i ~ c, then y^, . . . , ?/,;+fc„i are precisely the 
non-zero roots of the polynomial 

P^'^ = E P^Y' e C[Y] , 

0<j<n 

counted with multiplicities. 

The classical Newton-Puiseux theorem applies to the case where C is replaced by 
the field of (formal, or convergent) Puiseux series (a Puiseux series is of the form 
J2T=K a.kx^'' with afe G C, A' e Z and s £ N\ {0}), and shows |2=^n only In the 
classical statement of the theorem, the leading exponents Yi, are, up to an inversion 
and change of sign, the slopes of the Newton polygon, and the polynomials p^'' are 
defined in terms of the edges of the polygon. Since, when P = e{V), the graph of 
vexP is the symmetric, with respect to the main diagonal, of the Newton polygon, 
it follows from Proposition I2.1U| that the Yi and yi in Theorem 13. II coincide with 
the ones that are defined classically. 

Theorem 13.11 is a "precise large deviation" version of the Newton-Puiseux the- 
orem: we assume only the existence of asymptotic equivalents for the coefficients 
of 7^(6, ■), and derive the existence of asymptotic equivalents for the branches of 
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P(e, •). The Newton-Puiseux algorithm is sometimes presented for asymptotic ex- 
pansions, as in |Die68j . However, the equivalence between the two assertions of 
Theorem 13 . II does not seem to be classical. In particular, the asymptotics of some 
coefficients may be only known as being negligible: wc require that pi ^ only for 
those i such that (i. Pi) is an exposed point of the epigraph of P. 

Proof. We first provc[Tl=^|21 Let Q = (Y® Yi) • • ■ (Y©r„). Then, Q = Q, R((5) = 
{Yi<---< r„) and Q„_, = Fi ■ • • for alH = 1, . . . , n. Since yi{e), . . . ,y„(e) 
arc the roots of V{t, y) = counted with multiplicities, and Vn = 1, it follows that 
^(e,Y) =nr=i(Y-3^i(e))- Hence, {-iyv„-, is the sum of ah products 3^^, ■•■ J^^. , 
where are pairwise distinct elements of {l,...,n}. By the properties 

of "~" (stability by addition and multiplication) , and since 0^-^ j -^h ' ' ' ~ 
Yi ■ ■ 'Yi = Qn-i, we obtain that there exist po, . . . ,Pn~i G C such that Vj ~ pje^^ 
for all j = 0, . . . , 71 — 1. Putting p„ = 1, we also get 7^„ = 1 ~ Pn^'^" since Qn = 1. 
When i = l,...,7i — 1 is such that Yi < K^+i, • • • yi is the only leading term in 
the sum of all ■ ■ ■ yj., and then Pn-i = • • • y,; ^ 0. Moreover, for i = n, 

either F„ ^ (D, which implies that po = (— l)"yi ■ • ■ 2/ti 7^ 0, or y„ — (D, which implies 
that IVn = 0, Po = and Qo = 0. This shows that (ci, . . . , c„) = (Yi, . . . , y„) and 
P = Q are as in Point 13 

The remaining part of the theorem is obtained by a simple adaptation of the 
proof of the classical Ncwton-Puiscux theorem. When the Vj are only assumed to 
be continuous functions satisfying Point|3of the theorem, it follows from (|14I15I19|I . 
that e{V) > P, and since P >P= (Y ® ci) • • • (Y e c„), we get that e(P) > 
(Y©Ci) • • • (Y©c„). In addition, from H17|l and Point |2 of the theorem, we get that 

e(P)„_i = Pn-i = P„ci ■ C j for all i G {0, n} U {i e {1, . . . , n - 1} | c, < Q+i}, 

hence Lemma ICTl yields e(V) = P, therefore, R(e(P)) = R(P) = (ci < • • ■ < c„). 
Moreover, the first step of the Puiseux algorithm shows that, for all roots c 7^ (D of P 
with multiplicity k, there are exactly k continuous branches with leading exponent 
c. Indeed, when c = c.^ = • • • = Ci+^-i 7^ (D, the change of variable y = ze'^, and 
the division of P by transforms the equation 'P{e,y) = into an equation 

Q(e, z) = 0, where Q{-,z) extends continuously to with Q(0, z) = Since 
P(c) = PjC^ implies that n — i — k + 1 < j < n — i + 1, and since either i — 1 = 
or Ci-i < Ci, wc get that Pn-i+i 7^ 0, hence degp*-*' = n — i + 1. Similarly, wc have 
either i + fc — 1 = n or Ci^k-i < Ci+k- In the second case, we get pn-i-k+i 7^ 0, 
thus valp^') = n — i — k + 1. In the first case, i + fc— l=n, c = Cn, and po ^ 
or Po = (D. Since Po = (D implies c ~ Cn = (D, which contradicts our assumption, we 
must have po 7^ 0, hence again valp''^ = n — i — fc+l. Hence, degp^^'' — valp^'-* = k 
and the conclusion is obtained by the standard Lemma l3 . 21 below. Finally, if c = (D 
is a root with multiplicity fc, then valP = fc, Cn~k < Cn-k+i = ©, and pk =/= 0. This 
implies that (for all e > in a neighbourhood of 0) P(e, •) is a polynomial with 
valuation k, hence it has as a root with multiplicity k. □ 

Lemma 3.2. Let Q(e,Y) = X]"=o Qji^)^'' ' where the Qj are continuous functions 
of e G [0, eo) o,n-d let m = deg Q(0, •). Then, for any open ball B containing the roots 
of Q(0, •), there are m continuous branches Z\, . . . ,Zm defined in some interval 
[0, ei), with < ei < eo, such that Zi{t), . . . ,Z,„(e) are exactly the roots of Q(e, •) 
in B counted with multiplicities. Moreover, the roots of Q(e, •) that are outside B 
tend to infinity when e goes to 0. 
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Proof. We only sketch the proof, which is classicaL By the Cauchy index theorem, 
if 7 is any circle in C containing no roots of Q(e, •), the number of roots of Q{e, •) 
inside 7 is {2Tri)~^ dzQ{€, z){Q{e, z))~^ dz. By continuity of e 1-^ Q{e,-)j the 
number of roots of Q(e', •) inside 7 (counted with multiplicities) is constant for e' 
in some neighbourhood of e. Taking B as in the lemma, 7 = dB, and e = 0, we get 
exactly m roots of Q(e', •) in B for e' in some interval [0, ei). Consider now a ball 
Br D B oi radius R. For e' small enough, the number of roots of Q(e', •) in either 
Bji or B is equal to m, hence any root of Q(e', •) outside B must be outside Br. 
This shows that the roots of Q(e', •) that do not belong to B go to infinity, when 
e' — > 0. Finally, by taking small balls around each root of Q(e, •), with < e < ei, 
we see that the map which sends e to the unordered m-tuple of roots of Q(e, •) 
that belong to B, is continuous on [0, ei). By a selection theorem for unordered m- 
tuples depending continuously on a real parameter (see for instance |Kat95l Ch. II, 
Section 5, 2]), we derive the existence of the m continuous branches Zi, . . . , Z^. O 

Theorem 13.11 says that "the leading exponents of the roots are the min-plus 
roots" . 

Example 3.3. Consider V{e, Y) = + e^Y^ - e'^Y + e^^. The min-plus polynomial 
P ~ g{V) is the one of Figure^ hence its roots are ci = C2 = 3 and C3 = 7. We have 
— p(2) = Y'^ — Y and p'^''*^ = — Y + 1. Hence, V has 3 continuous branches around 
with first order asymptotics: yi ^ e^, ~ —e^ and 3^3 ~ e^. Theorem l3.1l statcs 
in particular that we need not know the asymptotic expansions of all the coefficients 
of r{e,Y): for instance, if P{e,Y) = + o{e^)Y^ - e^Y + e^^, the polynomials P 
and p^^\p^'^\p^'^^ are unchanged, so that we still have 3 continuous branches with 
the same asymptotics has above. 

Remark 3.4. If ^ G C"^" satisfies (|20f) . the characteristic polynomial of Ae: 
Vie, Y) ~ dct{YI—Ae) is an element of C[Y] , since C is a ring. Applying Theorem l3.1l 
to P, we can obtain, under some additional assumptions, first order asymptotics 
for the eigenvalues of A^. The difficulty is that the coefficients Vj of V need not 
have first order asymptotics (even if 7^ for all i,j) due to cancellations. Of 
course if the coefficients of Ae have Puiseux series expansions in e, the Vj also have 
Puiseux series expansions in e and a fortiori first order asymptotics. However, if we 
only assume that A £ satisfies (|20|) . we obtain that the Vj satisfy the con- 

ditions = 1 and Vj{e) ~ Pj^^' for some exponents Pj S Mmin computed using 
the exponents A^j (see Section IX^ . Hence, if the eigenvalues of A(^ have first order 
asymptotics. Theorem 13.11 gives the exponents of these asymptotics as a function 
of the Pj. 

3.3. Majorisation inequalities for roots of min-plus polynomials. The per- 
manent of a matrix with coefficients in an arbitrary semiring (5, O,®) can be 
defined as usual: 



where 6„ is the set of permutations of {1, . . . ,n}. In particular, for any matrix 



n 



per(A) = (g) 



cree„ 2=1 



■ 71 xn 
'mill ' 



71 
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and the formal characteristic polynomial of A is the polynomial 

n 

per(Y/e A) = (g)(Y,5,,(,) © e Mmi„[Y] , 

<T6e„ i=l 

where / is the identity matrix, and 5ij = 1 if i = j and dij ~ otherwise. The 
assoeiated min-plus polynomial function will be called the characteristic polynomial 
function of A. 

We next assume that A G C"^" satisfies (|20|l and that the eigenvalues CI [i = 
1, . . . ,n) of .4e have first order asymptotics, ~ A^e^'. We relate, in that case, 
the Ai with the roots of the characteristic polynomial of A. 

We need first to recall the classical definition of weak majorisation (see |M079j 
for more background). 

Definition 3.5. Let u,v ^ ^min- Let < • • • < ^(n) (resp. i'(i) < • • • < W(„)) 
denote the components of u (resp. v) in increasing order. We say that u is weakly 
(super) majorised by u, and we write u v, if the following conditions hold: 

"(1) ■ • • W(fc) > W(i) ■ • ■ V(k) Vfc = 1, . . . , n . 
In fact, the weak majorisation relation is only defined in |M079| for vectors of 
M". Here, it is convenient to define this notion for vectors of We also used 

the min-plus notations for homogeneity with the rest of the paper. The following 
lemma states a useful monotonicity property of the map which associates to a 
formal min-plus polynomial P its sequence of roots, R(P). 

Lemma 3.6. Let P,Qg Rmin[^] be two formal polynomial of degree n. Then, 

(22) P>Q and P„ = Q„ =^ R(P) R(0) . 

Proof From P > Q, we deduce P >Q. Let R(P) = (ci(P) < ■ • • < c,i(P)) and 
^{Q) = {ci{Q) < • • • < Cn{Q)) denote the sequence of roots of P and Q, respec- 
tively. JJsing P ^Q, Pn = Pn = Qn = Qn ^nd ®, WC get Ci (P) • • • Cfc (P) = 

P,i_fc(P«)"^ > Qn-kiQny^ = ci{Q) ■ ■ ■ Ck{Q), for aU k = l,...,n, that is 
R(P) R(Q). □ 

We shall also need the following notion of genericity. We will say that a property 
J^(y) depending on the variable y = (yi, . . . , j/n) G C" holds for generic values of 
y if the set of elements y e C" such that the property ^{y) is false is a proper 
algebraic variety. This means that there exists Q G C[Yi, . . . , Y„] \ {0} such that 
^{y) is false if Q{y) = 0. When the parameter y will be obvious, we shall simply 
say that ^ is generic or holds generically. It is clear that if and ^2 are both 
generic, then "^1 and ^2" is also generic. 

Since any polynomial q = I]n,...,i„6N *i,...,i„Yl' • • • Y,\" £ C[Yi,...,Y„] in n 
indeterminates can be seen as an element of C[Yi, . . . , Y„] whose coeflacients are 
constant with respect to e, we have: 

(23) e(<7) = Yi^ • • • Yj," G R,„in [Yi , . . . , Y„] . 

ii....,i„eN 
Q^^ .„#0 

We also define, for aU Y G M.'^i^: 

(24) qf.-*(Y) := Yi^ ■ • • Y;- G C[Yi, . . . , Y„] . 

ii,...,i„6N 
e(g)(yi,...,F„)=y;i---Y'^" 
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The following result is clear from the above definitions of e{q) and qy'', since 
when y / 0, ~ ye^ <;=^ y ^ ye^ ■ 

Lemma 3.7. Let q e C[Yi, . . . , Y„] and let Q = e{q) and qp"^ be defined by ^23\\ 
and 123}, respectively. Let y € C" , y e and Y e MJ,\i,j be such that ~ y^e^^ 
for i = 1, . . . , n. Then, 

(25) q{y,,...,y^)c,q^^'{y)e'^^''^ , 

and for any fixed Y , we have an equivalence in (|25|l for generic values of y ^ C" . 

Theorem 3.8. Let A £ C"'*" satisfy (|2U|I . Assume that the eigenvalues Cl, . . . 
of Ac (counted with multiplicities) have first order asymptotics, C\ ^ A^e^', and 
denote by A ^ (Ai < ••• < A„) the sequence of their exponents (counted with 
multiplicities) . Let T = (71 < • • • < 7n) be the sequence of roots of the min-plus 
characteristic polynomial of A. Then, 

(26) A r , 

and for generic values of a = (a^j) G C"'^", A = F. 

Proof. Since A — A^ £ Qnxn^ ^-^^ characteristic polynomial of A, Q(e,Y) := 
det(Y/ - Ae) belongs to C[Y]. Let Q = e(Q) S Kmin[Y] and let P = per( Y/ ® A) e 
I^min[Y] be the min-plus characteristic polynomial of A. By 114I15|I . e(Q) > 
per(YJ © e(-^)) and by (HOJ) and lO, e{-A) > A. It follows that Q = e(Q) > 
per(YJ© A) = P. Hence, from Lemma 1^ we get that R{Q) R(P) = T. More- 
over, by Theorem 13 . 1 1 applied to Q, we get that R(Q) = A, which finishes the proof 
of 

Let us show the genericity of the equality A = F. For all a E C"^", we consider 
the fc-th trace of a: 

tn{a)^ { Y sgn{a)Y[aj„(^jj 

Jc{i,...,n},#j=fe \creej jeJ 

For all A G M;\^„", we also set 

(27) tr,(^)= [0(g)^^^(^) 

Jc{i,...,n},#j=/c \<TGe,7jeJ 

Then, the coefficients of Q are given by Qfc(e) = (— 1)*^ tr„_fc(y^g), for fc = 0, . . . , 71 — 
1 and Qn = 1- The coefficients of P are given by Pk = tr„_fe(A), for fc = 0, . . . , ?i — 1 
and Pn ~ 1. By Lemma 13.71 we obtain that for any fixed (irreducible) matrix 
A G and any A G C"^" satisfying ^ with a G C"^" and A, trfe(A) ~ 

(trfc)^'**(a)e"''=('^^ for generic values of a G C"^". In particular, generically, Qfc(e) 
has first order asymptotics and e{Qk) = Pk, for all fc = 0, . . . , n. This implies that 
g = P, thus A = R(Q) = R(P) = F, generically. □ 

Remark 3.9. Since a result of Burkard and Butkovic |BB03j shows that we can 
compute the min-plus characteristic polynomial function of a matrix in polynomial 
time (by solving 0{n) assignment problems), Theorem 13.81 shows that the sequence 
A of generic exponents of the eigenvalues can be computed in polynomial time. See 
also |Mur9nilXBrT04) . 
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4. Critical values of min-plus matrices 

4.1. Schur complements and generalised circuit means. Wc now construct 
another sequence (3 = < • • • < /?„) using eigenvalues of min-plus matrices. 
First, we build by induction a finite sequence of min-plus square matrices and 
scalars S M, for 1 < ^ < A:, together with a partition Ci U • • ■ U C/c = {1, . . . , n}. 
We start with Ai = A. Then, for all £ > 1, we define 

(28) ai ^ Pn,in{At) 

and we take for Cg the set of critical nodes of Ag. We build, as long as CiU- • -UC^ ^ 
{1, . . . , n}, the min-plus Schur complement: 

Ai+i = Schur(C£, , y4f) . 

Due to the irrcducibility of A, Lemma 12.141 shows that Ag is irreducible, so that 
Ci ^ 0. Hence, the algorithm stops at some index k < n. By Proposition 12 . 1 5l 
we get that ai < • • • < a^. We call ai, . . . , the critical values of A. We define 
the multiplicity of the critical value ae as t^C^. Repeating each critical value with 
its multiplicity, we obtain a sequence (3 = (/3i < • • • < Pn) which will be called the 
sequence of critical values counted with multiplicities. 

Let us give now a graph interpretation of the exponents a/,. We set = and, 
for all £ = 1, . . . , fc, 

= Ci U...UCf, iV^ = {l,...,n}\C^-^ . 

For all paths p of G{A) and all ^ = 1, . . . , /c, we use the notations of Section 
and: 

IpIa :=bU - aiblci a£-iblcf-i , 

bl^ : = - IpIci \p\ce-i = \p\Ni ■ 

Proposition 4.1. The numbers ai defined in (|28fl satisfy: 

(29) ai = min -— j- , 

where the minimum is taken over all circuits c in G{A) which are not included in 
C^~^ . Moreover, c is a critical circuit of Ai if, and only if, there exists a circuit c' 
of G{A) such that c' n = c and c' minimises (|29|) . 

Proof. Using repetitively Lemma r2.14l we get that for all circuits c of G{Ai), \c\At = 
min |c'|^, where the minimum is taken over all circuits c' of G{A) such that c'ON^ = 
c. By the same arguments as in the proof of Proposition 12.151 we deduce the 
assertions of Proposition l4.1l □ 

Note that, as for Proposition 12.151 if c is included in C^^^, that is if the de- 
nominator in 129|) is zero, the numerator is necessarily nonnegative (by definition 
of ai-i). 

We say that a circuit c of G{A) is a critical circuit of order £ if |c|^ = af|c|^. 
We call critical graph of order i the graph G1{A) whose nodes and arcs belong to 
critical circuits of order £. Of course, G''{A) = G1{A). 

Proposition 4.2. We have 

(30) GK^)cG^+i(A) £^l,...,fc-l . 
which means that the nodes and arcs of Gi{A) belong to Gi^i{A). 
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Proof. If c is a critical circuit of order £, then by definition |c|^ = afjcj^. If in 
addition \c\^ 0, then cH TV*" = 0, whence \c\ce = and jcp^ = 0. It follows that 
I^Ia"^ = ~ C(e\c\ct = = Q;f+i|c|^+^, thus c is a critical circuit of order £ + 1. 
Otherwise, if |c|^ 0, then c minimises H29(> and since, by the arguments of the 
proof of Proposition 12.151 |cn N^Ia^ < |c|^, we obtain that c' = cO is a critical 
circuit of Ai. By definition of d, we get that the nodes of c' belong to Cg, thus the 
nodes of c belong to , which shows |c|^ = \c\ce or \c\^~^^ = 0. Since \cf^ ~ ae\c\^, 
we get |c|^^ = = |c|^+^, and c is a critical circuit of order £ + 1. □ 

Let £ G {1, . . . , fc} and let Di denote the min-plus diagonal matrix such that 
{De)jj = am if j £ Cm with m < £, and {Di)jj = ae if j E N^. For instance, if 
71 = 3, Ci = {!}, C'2 = {2,3}, ai = 2 and a2 = 4, then Di = diag(2,2,2) and 
D2 diag(2,4,4). We set 

Ai = D-^A . 

We also set 

G^(A) = G^.(A), i = ife, and7^ = i?fe. 

Lemma 4.3. 1/Fe have Ai = a£ Schur(C^~^, A^), for ( = 1, . . . ,k. 

Proof. We prove the lemma by induction on ^ = 1, . . . , fc. Since Ai = a^^ A and 
Ax = A, we get Ax = axM. If A^ = Schur(C'^~\ ^f), then using JTHll and (fTT|l . 

wc get 

Agj^x = Schur(C£, a^, Ai) = of^- Schur(C£, a^-^A^) 

= at Schur(Cf , Schur(C^-\ if )) = at Schur(C^ i^) 

= at Schur(C^ D tJ^xM+x) ■ 

Since {D'[^ De+x)jj ~ 1 for j G G^, and [DJ^ Dt+x)jj = a7^Q^£+i otherwise, it 
follows from ® that At+x = a^+i Schur(C^, Af+i). □ 

Proposition 4.4. For all I < £ < k, we have G'^{A) = C^lAi), At has min- 
plus eigenvalue 1, and the set of critical nodes of At is . Moreover, Gf {A) and 
G'^{A) n C^xC^ (that is the restriction of G'^{A) to the nodes of C^) have the same 
strongly connected components. In particular, G'^{A) = G'^{A) and all the nodes 
0/ {1, . . . , n} are critical for A. 

Proof. For all circuits c and for all £ = 1, . . . , fc, we get by Proposition 14. II |c|^ > 
at\c\^. Since, for aU circuits \c\j^^ = \c\a - ai|c|ci - • ■ • - at-x\c\ci^^ - at\c\Ne = 
\cYa ~ Q^^kl^j we get that p{At) > 0. Moreover, c is a critical circuit of order £ 
if, and only if, |c|^ = ai\cY , which is equivalent to |c|^^ = 0. This shows that 

p{At^ = and that c is a critical circuit of order £ if, and only if, c is a critical 
circuit of At- It follows that G^(A) = G'^{At). Since by Proposition l4.1l any critical 
circuit of At is of the form c' n where c' is a critical circuit of order t, the set 
Ci of nodes of G'^{At) is included in the set of nodes of G'i{A). Using (|3(J|I . we 
get by induction that C^ is included in the set of nodes of G'^(A). Conversely, 
since any critical circuit c' of order £ is such that c' n is a critical circuit of 
Ai, and since the set of critical nodes of At is Ci, the set of nodes of G'i{A) is 
included in ({1, . . . ,n} \ N^) U Ci = C^, hence is equal to C^. Finally it is clear 
that, by definition of At, G'^{Ai) C G"^(A), and since its set of nodes is C^, we get 
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G^iAi) C CiA) n C^xC^ . Conversely, since the restrictions oi A and Ae to C^xC^ 
are equal and since p{A() ~ p{A) = 0, any critical circuit of A with nodes in 
is critical for Ai. It follows that the strongly connected components of G^(A^) and 
G^iA) are equal. □ 

Example 4.5. To illustrate the computation of the critical values, consider 



(31) 



A 








OO 


oo 





oo 


1 


00 


1 


QO 


oo 


2 


oo 


00 


4 


5 



We have ai = 0, and the critical graph of A is composed of the circuit (1 
Thus, Ci = {1,2}. We have 

A2 = Schur(Ci,ai, A) 



!)• 



00 


2 ■ 




1 


00 




00 







00 


00 




■ 2 


2 " 


4 


5 


® 


00 


00 







00 




1 


00 




4 


5 



Hence, 02 — /3min(^2) ~ 2, with a unique associated critical circuit (3 — > 3), and 
C2 = {3}. (Recall our convention that Schur complements inherit their indices 
from the matrices from which they are defined, so that (^2)33 = 2 is the top left 
entry of A2.) We have A3 = Schur(C2, ^2, ^2) = 5 e 0*4 = 4, hence, = 4, with 
a unique associated critical circuit, 4^4, and C3 = {4}. 

To determine the critical graphs G'i{A), we use Proposition 14.41 which shows 



that Gj{A) = G%At). 
diag(0,0,2,2), and 



We already computed G\{A) = 6"= (A). Since D2 = 



A. 



D7,^A 



00 





00 


00 





00 


1 


00 


-1 


00 


00 





00 


00 


2 


3 



we deduce that G2{A) \ (A) consists of the circuit (1 
D4 = diag(0,0,2,4) and 



2^3-^1). Finally, 



(32) 



A3 = D-^A = 



00 





00 


00 





00 


1 


00 


-1 


00 


00 





00 


00 





1 



which shows that G'^{A) \ G2{A) consists of the circuit (3 
graphs are represented as follows 

1 2 

1 ^ ^ 2 ^ ^ 3 ^ ^4 



4^3). The critical 



Here, the graphs G'i{A), for £ = 1, 2, 3 are represented in black, magenta (medium 
gray), and green (light grey), respectively; for readability, a node or arc is drawn 
with the colour of the minimal graph G|(A) to which it belongs. 

For such a small example, the critical circuits could be obtained by mere inspec- 
tion. In general, G|(A) = G''{Ai) can be computed in polynomial time thanks to 
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Proposition which shows that G'^{Ai) coincides with the union of the strongly 
connected components of Sat{A£, V), for any eigenvector V of . 

4.2. Majorisation inequalities for critical values. We now state a second ma- 
jorisation result, which should be compared with Theorem 13.81 

Theorem 4.6. Consider an irreducible matrix A £ R^^^- Let F = (71 < • • • < 7„) 

be the sequence of roots of the min-plus characteristic polynomial of A and let (3 = 
(/3i < • • • < I3n) be the sequence of critical values of A, repeated with multiplicities. 
Then, 

(33) T 13 . 

Proof. Let P = pcr(Y/ © A) be the min-plus characteristic polynomial oi A,T = 
R(P) and Q = (Y ® /3i) ■•■(¥© Let V be an eigenvector of A (for instance 
any column A* , since by Proposition ^31 p{^) = 1 ^^^d all the nodes of {1, ... , n} 
are critical). Let W — diagV. Since AV — V, we get W~^AW1 = 1, where 
1 is the vector with all entries equal to 1. Therefore, W~^AW1 = Dl, thus 
{W-'^AW)^j > for all j = 1, . . . , n. Using we get 

(34) trfc(A) = UuiW'^AW) >pi---f3k ■ 
Then, 

P = Y" © tri(A)Y"~i ® • ■ • © tr„(A) > Y" © /3iY"^i © . . . © /^^ . . . /3„Y° 

= (Y©/3i)-..(Y©/3„) = Q . 

From LemmaESl we deduce R(P) R(Q) and since F = R(P) and R(Q) = /3, we 
obtain ((23 • □ 

Wc next characterise the cases where the equality holds in (j33|l . We say that a 
graph G has a disjoint circuit cover if there is a disjoint union of circuits containing 
all the nodes of G. This property, which is equivalent to the adjacency matrix of 
G having full term rank |BR91I Section 1.2], can be easily checked: it reduces to 
find a perfect matching (or to compute a matching of maximal cardinality) in a 
bipartite graph. 

Theorem 4.7. Consider an irreducible matrix A £ Rmin'^- F = (71 < • • • < 7„) 

be the sequence of roots of the min-plus characteristic polynomial of A, and let (3 = 
(/3i < • • ■ £ /3ri) be the sequence of critical values of A repeated with multiplicities. 
For all i 6 {1, . . . , A:}, where k is the number of critical values of A, the following 
assertions are equivalent: 

(1) 7, = forj e + 1, . . .,#G'}, and 71 ■ ■ •7#c*-i = Pi ■ • ■/3#c*-i; 

(2) G£_j^(j4) and G'^{A) have a disjoint circuit cover. 

In Theorem 14. 71 we use the convention that Gq{A) is the empty graph and that 
it has a disjoint circuit cover. Recall also that G° = 0. 
The proof of Theorem 14 . 71 relies on the following lemma. 

Lemma 4.8. The equality 

(35) tT#c4A)=Pi---f3#c^ 
holds if, and only if, G^(j4) has a disjoint circuit cover. 
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Proof. Let us first assume that G|(A) has a disjoint circuit cover. Since by Propo- 
sition ^31 the set of nodes of G|(^) is , there exists disjoint elementary circuits 
ci, . . . ,Cq in Gg{A) which cover all the nodes of C^. Let a be the permutation of 
the nodes of which consists of the circuits ci, . . . , Cq. We obtain, using H27|) : 

tr^c^A) < (g) = /3i • • ■ /3#c' (^f)j-O) = A ' ' ' P#c' , 

since, by Proposition ^31 ci, . . . , Cq are critical circuits of Ag and piAi) = 1. Since 
it follows from H34|l that tr^(je{A) > (3i ■ ■ ■ P^ch we have proved (|35|1 . 

Conversely, let us assume that (|35f) holds. Let W be as in the proof of Theo- 
rem 14.61 By (|27|) . there exists disjoint circuits ci, . . . ,Cg of G{A) such that |ci| + 
••• + |c,| =#C^ andtr#c^(A) =0^.g^^^...^^^A,,(,) ={g)^.g^^u...u^^(M^-iAM^),,(^ 
where a is the permutation of the nodes of consisting of the circuits ci,. . . , 
Cq. Since W~'^AWl > Dl, we obtain that iT^c^iA) > 0jeciu - uc, ^ii- If 
ci U • • • U Cg 7^ C*", we obtain, using /3„ > • • • > /3^c'^+i > /3#c^ > • • • > 
that tr-^pf (^) > /3i • ■ • P^c'^i ^ contradiction. Therefore, ci U ■ • ■ U = C^, and 
since /3i • ■ • /3#c'' = tr#c''(^) = <2)jec'' ^j't(j) , we get 0j-6ciU-uc, = 1- 
Since = 1, the circuits ci, . . . , Cg, which are critical for Ai, are critical circuits 

of G'i{A) (by Proposition 14. 4|) . Hence, G^(A) has a disjoint circuit cover. □ 

Proof of Theorem \4-7\ Let P = per(Y/0 A) be the min-plus characteristic polyno- 
mial of A. By Lemma [2 .111 we have 

(36) Pn-i = 7i ■ ■ ■ 7i < Pn-i = tri(A) , with equality when ji < 7^+1 . 

We prove 13 =>n] Assume that Gj_i{A) and G^{A) have a disjoint circuit cover. 
Combining the inequality in (|36|) with H35|) . we get 71 • ■ ■ 7;^c* /3i ■ ■ ■ P^c"- Sim- 
ilarly, 71 ■ • •7^(7«-i < /3i • • ■ P^c'^-^- Using (|33|l . we get the reverse inequalities 

(37) 71 ■ ■ ■ 7i > /3i ■ • ■ P], for j = 1, . . . , n 
so that 

(38) 71 ■ • ■ 7#cf = /3i • • ■ /3#c* . 

(39) 71 • • •7#cf-i = /3i • • •/^#c'!-i ■ 
Dividing ijSHl) by (1251), we get 

(40) 7#cf-i+i ■ • • l#c' = /3#c*-i-i-i • • • /3#c* = af^' 

(recall that #C£ = - #C^'^). Taking j = #G^-i 1 in JSTjl, and using 
we get 

7#c«-i+i > /3#c«-i+i = «£ • 

Since (7^) is nondecreasing, 7^ > ^^c'^-^+i — holds for all j £ {^C^~^ + 
1, . . . , #G^}, hence, if 7^ > for some j S {#C^"^ -I- 1, ... , #G^}, wc would 
have 7#c*-i+i " ' ■7#c* > oif^\ contradicting (^UJ). Therefore, 7^(74-1+1 = 
7#c* = = /3#cf-i+i =:••■ = /3#c*- 

We next prove[T]=^El By assumption, (jSHJ and ^ hold. Taking j = + 1 
in l|57|l and using we have 7^c'«4.i > /3^c''-i-i- Since P^^i^i > f3^c' ~ l#ch 
we have j^c'+i > 7#c' i so the equality case in H36() yields 

(41) 71 ■ • •7#c* ti-#c«(^) • 
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Taking now j = ^ — 1 in (pITji . and using we get P^c'-^ ^ 7#c*-i i hence, 

T^c'^-i+i = /3#c'-i+i > /3#c*-i 7#c*-i: s-nd the equahty case in yields 

(42) 71 ■ • • 7#c*-i = tr^c-'-i (A) . 

It follows from Lemma Ol and from (j^HJ, EB) and that G'^(A) and 

G^_i(A) have disjoint circuits covers. □ 

Corollary 4.9. If G'i_-^{A) and G'i{A) have a disjoint circuit cover, then ai is a 
root of multiplicity ^Ci of the min-plus characteristic polynomial of A. 

Proof. Since 7j = [ij = for j E [jjzC^^^ + 1, . . . , #C^}, ai is a root of multiphc- 
ity at least — #C^^^ = of the characteristic polynomial of A. Moreover, 
we showed in the proof of =^ 0' of Theorem 14.71 that 7^(7?+ 1 > ^^c' and 
7#C''~i+i > "f^c'-'^- Thus, is a root of multiplicity exactly of the charac- 
teristic polynomial of A. □ 

5. ASYMPTOTICS OF EIGENVALUES 

5.1. Statement and illustration of the result. We next show that under some 
non-degeneracy conditions, the first order asymptotics of the eigenvalues of Ae are 
given by the critical values of A. If G is any graph with set of nodes 1, . . . , n, and 
if & G C"^", the matrix 6*^ is defined by 

I otherwise. 

Let G be cither the critical graph of A or the saturation graph Sat (A, V), for any 
eigenvector V oi A (since by Proposition 14.41 all the nodes 1, . . . ,n belong to the 
critical graph of A, we can take for V any column of A*). 

We construct the following conventional Schur complements: 

(43) 5^=0'^, / = Schur(C^-\si), £ = 2,...,A: . 
The Schur complement s^ is well defined as soon as the matrix 

(44) = agf_i 

is invertible (we adopt the convention that r^ is the empty matrix, and is invertible). 
We shall also need the following matrix: 

(45) t' = • 

When both and s^^^ are well defined, t^^^ is invertible and we can compute 
from s^^^ thanks to Hll() : 

s*" = Schur(C£_i,s^-i) . 

We say that a function of e, /(e), is of order w(e") if limg^o !/(£)£ "I = +oo- 

Theorem 5.1 (Generalised Lidskii-Visik-Ljustcrnik theorem). Let , r^, t^, £ = 
1, . . . , /c be constructed as in H43I44I45() with G = G'^{A) or equivalently with G = 
Sat(A, V) for some eigenvector V of A. Assume that the matrix r^ is invertible for 
some I < £ < k, and let X{, . . . , denote the non-zero eigenvalues oft^ (here and 
in the sequel, eigenvalues are repeated with multiplicities). Then, the eigenvalues 
of Ae can be grouped in 
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(1) mg eigenvalues with asymptotic expansions 
(46) /:f■^~A^e"^ l<j<m, 



(2) ^ eigenvalues of order ui{e'^^), 

(3) #-/V^ — mi eigenvalues of order o(e"*). 

In particular, when t^, . . . jt*^ all are invertible, for all 1 < £ < k, Ae has exactly 
eigenvalues of order e°'' , whose asymptotics are given by (|46() . 

Wc prove Theorem 15. II in Section 

By Proposition the saturation graph Sat(A, (defined in Section and 
the critical graph G'^{A) have the same strongly connected components. This ex- 
plains why, in Theorem l5.1l one can use either the graph G = G'^{A) or the graph 
G = Sat(i,F). 

The following result, that we also prove in Section shows that the assump- 
tions of the theorem are generically satisfied, if we assume that the critical graphs 
have disjoint circuit covers: 

Proposition 5.2. Let I ^ l,...,fc. Assume that G'i_-^{A) and G1{A) have dis- 
joint circuit covers. Then, r^ and t^ are generically invertible, so that the number 
of eigenvalues of A^ with an equivalent of the form Ae"*, where A G C \ {0}, is 
generically =f^Ci. 

Example 5.3. To illustrate Theorem 15. II consider the matrix 



(47) 

so that {Ae)ij 



1 

e 

^2 



e 

^-2 



' , with 





' 1 


1 


1 " 




1 





4 


a = 


1 


1 


1 


, A = 


oo 


1 


-2 




1 


1 


1 




1 


2 


oo 



We have Pmin{A) = —1/3, and G'^{A) consists of the critical circuit: 

-2 
1"-^ ^2 ^ 3 



so that the construction of the critical classes stops with Ci = {1, 2, 3} and k = 1. 
Then, G^^^A) = G^A) = G^iA) covers all the nodes (see Proposition 14. 4|l. hence, 

1 

1 

1 

Since the spectrum of s^ is j^}. Theorem 15 . II shows that the spectrum of Ae 
consists of the three eigenvalues 



-1/3 r3 
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Example 5.4. To give an example in which different exponents appear, consider 

ai2 



A, 



021 



£023 



6^034 

e^a44 



where G C, and "•" denotes a zero entry. The associated matrix of exponents 
A is given by 1)3 1() , and we saw in Example 14.51 that the critical values of A are 
ai = 0, a2 = 2, as = 4, with Ci = {1, 2}, C2 = {3}, C3 = {4}. The critical graph 
G — G'^{A) of the matrix A = A3 of 1)32(1 was represented in Example 14.51 Thus, 



-.G 



ai2 



0-21 

031 



023 



043 



034 



The eigenvalues of the matrix = "q^] are the square roots of 012O21. Let 
us assume that 012O21 ^ 0. Then, Theorem 15 . II shows that Ae has two eigenvalues 
with asymptotics of the form ^ ^, where = 012O21. Moreover, 



= Schur({l,2}, 



O31O2I "^23 
O43 



034 



.2 _ -1 

t — —a3ia2i 023 



If we assume additionally that 031023 7^ 0, Theorem 15.11 shows that Ae has an 
eigenvalue with asymptotics —0310^]^ 0236^. Finally, as soon as the matrix 
is invertible, i.e., as soon as detr'^ = 012023031 7^ 0, the Schur complement 



043023^ 02103]^^^ 034 



= s^ = Schur({l,2,3},si) 

is well defined. (When is well defined, that is 012021 0, and 031023 7^ 0, 
we may obtain equivalently as Schur({3}, s^).) Thus, when 012O23O31 ^ and 
043O21O34 7^ 0; Theorem 15.11 shows that Ae has an eigenvalue with asymptotics 

r -1-14 
i-e ~ 043023 02l03j^ 0346 . 

5.2. Proof of Theorem lS.ll and Proposition l5.2l For the proof of Theorcm l5.1l 

we need to use the following lemma, which follows readily from the definition of 
determinants. 

Lemma 5.5. Ifb,b G C"^" have two digraphs G{b) and G{b) whose circuits (or 
equivalently, whose strongly connected components) are the same, and if bij = bij 
for all arcs {i,j) belonging to circuits of G{b) or G{b), then, det 6 ~ det&. 

Let V be an eigenvector of A and let Sat = Sat(A, V). The change of variables 
A = fie"'' , for some 1 < i < k, transforms the characteristic polynomial of Ae into 

dctifie"' I ~ Ae) = detie^") det{ne"'e^^' I ~ e^^' Ae) = dct{e^')V{e, fi) 

where P(e, /i) ^ dctific'^'e"^' I ~ e°^\'^'^s(vr' j,^^di^s{V)-j 

If C C £ arc finite sets, we denote by E^, the L x L diagonal matrix such that 



{Eh)^ 



1 for i G C , 
for i G L \ C . 
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If L = {1, . . . , 7i}, we shall simply write Ec instead of i?^. We have 

(48) e'^'e'^^' ^Ej^i , 

hence ■p(e, /i) — >e^o ''^(O, m), where 

7^(0, /i) = det(AiSjvf - £;c*aS''*) . 

Since Sat and G^(A) have the same strongly connected components (by Proposi- 
tion I^J, Lemma [5.51 vields : 

P{0,ti) = dct(^£;jvf - Ecea'^"'^^^) . 

The same arguments also show that the invertibility of the matrix is independent 
of the choice of G = Sat or G = G'^{A) in H43|l . Hence, if s^, r^, are constructed 
as in (|43I44I45|) with either G = Sat or G = G^iA), and if is invcrtible, then 

7^(0, m) = M*"^'^' detifiESl - ag^.c*) = M^"^'^' det(-/) det{fil - t') . 

From Lemma [3.21 applied to V{e,\), there exists #iV^ continuous functions e > 
£™'-', with j = 1, . . . , #Cm and m = £, . . . ,k, such that are the roots of 

7^(e, Y) for all e small enough. Hence, Cq-' are the eigenvalues of and = 
for m > £. The other roots of 7^(6, Y) tend to infinity. This shows Theorem 15. II 

We finally prove Proposition 15. 21 If the set of nodes of G^_]^(^) can be covered 
by disjoint circuits, it follows from Proposition 14.41 that these circuits also belong 
to G'^{A) n G^~^ X C^~^. By definition of , for generic values of a = (fly ), these 
circuits belong to the graph of r^, which implies that the determinant of is 
gcnerically non-zero. Thus, is generically invcrtible. The same argument shows 
that if GiiA) can be covered by disjoint circuits, r^"^^ is generically invertible, and 
since = Schur(G^~^, r''"'"^) is the Schur complement of the generically invertible 
C^~^ X C^~^ submatrix of r^+^, namely r^, in the generically invertible matrix r^"*"^, 

must also be gcnerically invcrtible. Thus, — #Gf generically in Thcorem l5.1l 



6. ASYMPTOTICS OF EIGENVECTORS 

6.1. Statement and illustration of the result. We now consider eigenvectors. 

Theorem 6.1. Let s^, , t^, £ = l,...,k be constructed as in Theorem \5.1[ 
Assume that the matrix is invertible, for some 1 < £ < k, that /j. ^ is a simple 
eigenvalue of t^ , and let V be any eigenvector of Ai. Then, the equation 

(49) (^i^ATf -aS"*(^^'^))w = 0, 

has a unique solution w = {wj) G C"\{0} up to a multiplicative constant. Moreover, 
there is a unique eigenvalue with asymptotics ~ M^"' ; '"^'^ if Wi ^ 0, any 
eigenvector Ve associated to this satisfies {ye)i ^ for e small enough, and 

(50) , forjG {!,...,«} . 
We prove Theorem 16. II in Section IfT^ 
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Example 6.2. To illustrate Theorem l6.ll let us pursue the analysis of Example 15. 31 
We already showed that the eigenvalues of the matrix H47|l have asymptotic equiv- 
alents of the form ^e^^/^, where ^ is a cubic root of 1. When fi ~ ^, any solution 
of H49|) (with £ = 1), is proportional to w = [1, ^, ^^]"^. Since A has a unique critical 
class, Ci = {1,2,3}, by Theorem 12.31 A has a unique eigenvector, up to a scalar 
factor, and we can take V = [0,-1/3,4/3]"^ = A*-^^. Theorem 16. II shows that any 
eigenvector Ve associated to the eigenvalue ^e^^/'' is equivalent to 

[l,e.-V3_^2^4/3]T 

up to a scalar factor. 

When Wj = 0, Theorem 16.11 gives a poor information on the asymptotics of 
{Ve)j. Moreover, when A^ has several critical classes (so that the eigenvector V is 
non unique) the non-zero character of Wj depends in a critical way of the eigenvector 
V which is selected. 

Example 6.3. The following example illustrates the importance of the choice of the 
eigenvector V in Theorem 16. II Consider 



1 

-2e 



2e^ 



which is such that {Ae)ij — aije^^' with 



" 


1 


3 




1 


1 


1 ' 


1 


2 


oo 


, a = 


-2 


1 




3 


oo 


2 




1 




2 



We have ai 
and 



A = 



Pmin{A) = 0, with a unique critical circuit (1 — > 1). Hence, Ci = {1}, 



A2 ^ Schur({l},yl) 
Thus, a2 = 2, C2 = {2, 3}. We have 



A 







-1 



1 


00 



3 

00 




Since the critical graph of A, which is the union of the complete graph on {1,2}, 
and of the loop (3 3), has two strongly connected components, {1,2}, and {3}, 
the eigenspace of A is spanned by the two vectors A*^, i = 1,3. Let us take 



V = A*3 = [3,2,0]^ , 

for which the saturation graph is obtained by adding the arc (1 
graph of A. Taking G = Sat(A, V) in (gS)), we get 

1 1 1 



3) to the critical 



(51) 



a 



-2 1 



Since = 1, Theorem 15.11 shows that has a root with asymptotics /3c ~ 1, and 
since = = Schur({l}, s^) = [j^ ^] has roots 2,3, Ae has two eigenvalues with 
respective asymptotics Ce ~ 2e^, and Ce ^ 3e^. Let us compute for instance the 
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asymptotics of the eigenvector Vc associated to Lt ^ 2e^, using Theorem 16 . II with 
1—2 and fi — 2 (thus A2 ~ A). With the previous choice of V, we need to solve 
the system (I49II . which, by H51|l . speciahses to 

wi + W2 + W3 = 0, —2wi — W2 ~ 0, = . 

All the solutions of this system are proportional to w ~ [1, ~2, 1]"^. Thus. Theo- 
rem shows that up to a multiplicative constant, 

Consider now the alternative choice of V: 

V = i*i = [0,-1,1]'^ . 

Then, Sat(42, V) is obtained by adding the arc (3 — > 1) to the critical graph of A. 
Theorem 16. II yields that {Ve)i — Wie^', where 

Wi + 102 — 0, — 2wi ~ W2 — 0, Wi = 

and since all the solutions w are proportional to [0, 0, 1]"^, we learn only from H5U|) 
that (Ve)i/(Ve)3 ~ Oe^^, and (Ve)2/(Ve)3 ~ Oe"^, a very poor information. 

Remark 6.4. When n is not a simple root of t^, the first order asymptotics of the 
eigenvector may be ruled by higher order terms in the expansions of the entries of 
Ae, see |ABG98| for a special case. 

6.2. Proof of Theorem 16.11 We first observe that by Theorem 15. II there is only 
one eigenvalue of Ae equivalent to fj.e°"^. Then the associated eigenvector, Ve, 
is unique, up to a multiplicative constant, since for e small enough, is a simple 
eigenvalue of A^. 

To prove Theorem 16.11 we perform the change of variables Vc = e'^^^^^We and 
Ce ~ Me^"', where AA^ — * 1^ when e ^ 0. After multiplying Vc by a constant, we 
may assume that J2i<j<n ~ 1- From A^Ve = C^Ve, we get 

where eC7igdiag(v)-i_4^gdiagy _^ QSat(^,y) ^^^^ ^ _^ q Together with this 
implies that any limit point w of Wc when e ^ satisfies 

(52) a^''^^^''^^w ^ fiENcw, and |uii| + •• • + = 1 . 

To show that the solution w of (|52|l is unique, up to the multiplication by a complex 
number of modulus 1, we shall prove that ^Ej^e — aSat(Af,v) ]^,^ rank n — 1. 

Since, by Proposition 12.41 Sa,t{Ae,V) and G^iAi) = G^(A) have the same 
strongly connected components, applying Lemma [5.51 to the matrices b = b{X) = 
XEj^j - aSa*(^*'^) and b = 5(A) = \Ej^, - a<=?(^), with A G C, we get det6(A) = 
det6(A). Moreover, since G'^{A) and the restriction of G'^(A) to have the same 
strongly connected components (see Proposition I4.4|l . then by Lemma 15.51 asrain. 
det6(A) = det(A£;g^'' - r'^+i)A#^'+' = det(-r'') det(A/ - t'^)X*^'^\ which yields: 

(53) det5(A) = det(-/)det(A/-t'')A*^'^' ■ 

Hence, detb{^) ~ since /i is an eigenvalue of t^, and fiEj^e ~ Q,Sat(A{,v) -^^^^ rank 
< n. Since is a simple eigenvalue of and ^ 7^ 0, is a simple root of the 
equation det 6(A) = 0. Hence, the partial derivative d\ det 6(A), evaluated at A = ^, 
is non-zero, which implies that there is a subset L of {1, . . . , n}, of cardinality n — 1, 
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such that det(6(/i)i.i) ^ 0, which shows that [lEj^i — aSat(yi^,v) j^g^g rank n — 1. 
Thus, (|49|l has only one non-zero sohition, up to a scalar multiple, which implies 
that all the solutions of 152|) are of the form Cw, where C g C is such that |C| = 1, 
and w is any solution of 152|) . Let us pick i such that iti^ ^ 0. Since all the limit 
points of We are of the form C,w, with \C, \ = 1, we get (yVc)j/(yVc)i Wj/wi when 
e ^ 0, and since e'^'^'s^We, we get (jSOJ. 

6.3. On the choice of the eigenvector V . We now show that there is, in some 
sense, a canonical choice of V in Theorem 16.11 Denote by Cf, . . . ,C5j the critical 
classes of Ag^ and by C|,...,C£* their restrictions to Ci. By Proposition 14.41 
cf , . . . , C^^ are the strongly connected components of G"^(4)nC^ x and they cover 
. Moreover, one can deduce from Proposition 14. II that ior ~ 1, i/i, Cg is 
either the empty set or a critical class of the matrix Ai , and that U • • • U C^' ~ Cg . 
Then, when is invertible, the characteristic polynomial of can be factored as 

(54) det(A/ ~t^)^ Q\ (A) • • • Q^^ (A) 

where Q^(A) = det(A/ - t^c",,C",) if ^ ^ and Q^(A) = 1 otherwise. Indeed, 
taking G — G'^{A) in (|43|) . using the fact that Cf , . . . , C^^ are the strongly connected 
components of G'^{A)r\C^ x C^, and using the block triangular structure of XE^^i — 
aG'^(A)^ we get 

(55) det(-/) det(A/ - = det(A£;g^'' - a%^^}) 

= ndct(Ai?g-a5(,f) 

= det(-rO n det(A/-<S,^cr) • 

l/=l,...,l/f , C^^7^0 

Since is invertible, this shows H54|) . Thus, if /.( 7^ is a simple root of det(A/ — i''), 
there is a unique v e {1, . . . , i/'^} such that /i is a root of the polynomial (5^(A). 
Denote by i/(/i) this index. Let be an eigenvector of A^, for instance V = 
(A£)*j with i G C^. By the same arguments as in the proof of (|^ . one can 
show that (|55|l remains valid if we replace C^lAf) by Sa,t{At,,V). Hence, for any 
V ^ {^iEpfi — a^^'^^^''^^)ci .ci is invertible. Moreover, since ^ 7^ 0, {pLE^^i — 

QSat(Af,v')-j^^^^^^^^ invertible. One can then deduce, using the block triangular 
structure of ^Ep^t _Q,Sat(A;,v)^ ^^i^^t if there is no path from i to ^^(^j in Sat(A^, V), 
then Wi = 0. In particular, using Proposition 12.61 one deduce that MV = (^£)*j 
with j ^ \ C'^(^), then there exists a final class of Sa,t{Ai, V) different from 
hence Wi = for all i ^ Cf,. This observation explains Example 16.31 and 
it also suggests that the choice V = i^e)* j with j G ^^(^^ is canonical (note that 
different choices of j G C^i^^-j yield proportional vectors V). However, in the case 
of eigenvectors, there does not seem to be a simple analogue of Proposition 15.21 
(characterising the cases where generically w has non-zero entries). 
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7. The theorem of Visik, Ljusternik, and Lidskh revisited 



7.1. statement of the theorem. We now show that the theorem of Visik and 
Ljusternik |VL60| and Lidskii |Lid65j can be obtained as a coroUary of Theorem l5.1l 
and that Theorem 15. II allows to solve cases to which the classical result does not 
apply. The presentation of this subsection is inspired by |MB()97j . that the reader 
may consult for a general discussion of the theory of Visik, Ljusternik, and Lidskii. 

Lidskii )Lid65j considers a matrix of the form Ae = Ao + eb, where b G C"^" and 
Ao G C"^" is a nilpotcnt matrix. Wc shall need specific notations for Jordan matri- 
ces. Let N[q] denote the <Z x q nilpotcnt matrix such that {N[q])ij = 1 if j = i + 1, 



and iN[q]), 



otherwise. For m > 1, we define — N{q)- 



-Niq) 



(m-times), where + denotes the block diagonal sum, and, given a decreasing se- 
quence gi > ^2 >•••>%> 1, and mi,...,TOfe > 1, we define iV[™j;;;;;g"'' ] = 
^[Z'] + ---+N["qk]- For instance, when qi = 3,toi = 1,^2 = 2,m2 = 2,^3 = 
1 , 7713 = 1 J we have 



(56) 



' ■ 1 ■ 
■ ■ 1 
□ • • 


□ • 


□ • 


□ 


□ • • 


• 1 
□ ■ 


□ • 


□ 


□ • • 


□ ■ 


• 1 

□ • 


□ 


_ □ • • 


□ ■ 


□ • 


□ 



where, again, "•" represents (why some zero entries are written inside boxes will 
be explained below). We consider the case where is equal to ^^["i';:::;^" ] . If 
1 < i < k, we finally define the {nii -!-••■ + 77i£) x (7771 + • • • + 777^) submatrix 
of b, obtained by considering only the bottom rows and first columns of the Jordan 
cells of sizes qi x qj, i,j = l,...,£. For instance, in the case of H56(l . 



$1 = [ 631 ] 



bji 



674 



^36 
^56 



*3 



bsi 
b5i 
bji 
bsi 



b34 
674 



^36 
^56 
bjG 

bse 



bas 
658 
678 

&78 



The corresponding positions in the matrix ^0 were depicted by boxes in (|56() . 
convention, $0 is the empty matrix, and is invertible. 



By 



Corollary 7.1 ( jLid65l Th. 1]). A ssume that both <i>£„i and are invertible, for 
some 1 < i < k, and let Ai,...,Amf denote the eigenvalues 0/ Schur($£_i, $1). 
Then, Ae has miqi eigenvalues with asymptotics 

Ce £,e'^^'^' , where = A, and i = 1, . . . , 777^ 

(for each Ai, all the qe-th roots ^ of Xi are taken). 

Of course. Corollary 17.11 can be stated in an equivalent "coordinate free" way, 
by using left and right eigenvectors associated to the different Jordan blocks, 
see |Lid65j . In fact, Moro, Burke, and Overton observed that we need not re- 
quire to be invertible in Corollarv l7.1l when is sing ular, IMlj()97l Th. 2.1] 
shows that to each eigenvalue Ai G C of Schur($^_i, $£) corresponds q^ eigenvalues 
of Ae with asymptotics jC^ = ^e^^'' + o(e^/''«) where = A^. 
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7.2. Derivation of Corollary 17. IL Let us denote by ^[ ] the matrix 

of exponents associated to = N ["ql]',',',]^^^ ] +eb: A[' 

'Iju'.'.'.l^k'' ] is obtained from 

N [ "ql'.'.'.'.fi^^.'' ] by exchanging zeros and ones. For instance, ^ [ 2 ] = [ i i ] corresponds 
to — [q q] + e6. The following lemma is straightforward. 



Lemma 7.2. Let qi > • • • > (/fc > 1, irii, . . . ,mk > 1. The matrix A\^"qll',',',[™^'' '\ has 
min-plus eigenvalue 1/qi, and set of critical nodes Ci ~ {1, . . . ,7711(71}. Moreover, 

(57) Schur(Ci,iM,^["^i;::::r]) = 



It follows from Lemma |7.2I and in particular, from the recursive property (|57|l . 
that the sequence of critical values of ] is (ai, . . . , 0;^) — {1/qi, . . . ,1/qk), 

and that the associated critical classes are Ci = {1, 7771(71}, .... C\ — 
{^i^i'mgqg + 1, . . . ,'^g^^m^qi}. Recall that the diagonal matrix D is defined 
from the ai and Ci. 

An eigenvector T^[™J^;;;:9"'° ] of D^^ A["qll',','/J,^^''] can be built as follows. For all 
q > 1, we set V[q] = [0, 1/(7, . . . , ((7 — l)/q]'^, then, for 771 > 1, we define ] = 



V 



-V 



(777-times), where 



denotes the concatenation of vectors, and, 
finally, we set F[™J;;;;:g"'° ] = "/] + ■•• ] • I* is easy to see that V = 
V[^^l]'.'.'.'J^^'' ] is an eigenvector of A = ] , and that the corresponding 

saturation graph is the union of the graph of ^^["gi"!!;^^* ] and of the arcs {i,j), 
where i is the index of a bottom row of a Jordan block of N ["^{[[[[[ql'' ] , and j is 
the index of the left column of a Jordan block of A^[™i:!!!^g"'" ] . Since Sat(A, V) is 
strongly connected, it is also equal to G'^(A). For instance, for Ac = A^[3'2'i] + ^6, 
and G = G%A) = Sat(i, V), we get 



(58) 



• 1 • 
■ 1 

631 . . 


634 • 


636 • 


638 


651 • 


• 1 

654 • 


656 • 


658 


671 ■ ■ 


674 • 


1 

676 • 


678 


&81 ■ ■ 


^84 • 


686 • 


688 



The statement of Corollary 17.11 becomes a special case of the statement of Theo- 
rem provided the following identity is proved: 



det(A/ - t^) = det(A« - Schur($(?_i, $f)). 



1, 



This can be seen immediately by noting that t^ is a matrix of cyclicity , which can 
be put, by applying a transformation 1— > P^t^P^^, for some permutation matrix 
Pi, in block circular form 



(59) 



Pit'p; 



'mf(gf-l) 



Schur($£_i, 

where Iq is the identity matrix of order q, and where the "•" represent blocks with 
values. 

Indeed, by and we get: 



(60) 



Schur(C^ 



-1 ^ 
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and for each £ = 1, . . . , fc, there exists a matrix corresponding to a permutation 
of preserving Ci, such that in block form we get: 



1 _ 




■■■+me-i{qe-i-l) 








$22 



where $/> 



and — $£-1 (for each £, the indices of correspond to 

the nodes of Ce = {Y^tZi ^iQi + 1, • ■ • , J2i=i "i^gi} of the form X^ti + 
with m = 1, . . . , mf). Hence, taking for Pi the restriction of Qt to C^, and using 

the fact that *] ^ = ^-i * ^ for all invcrtible matrices and $, we get (|59|l . 

For instance, in the special case of H58|l . and ^ ~ 2, we get 



Q20'C2 C2Q2 ^ — 





1 


1 












634 


636 • • 


1 _ 








1 • 
• 1 




651 ■ 




654 


656 • • 




. 671 ■ 




674 


676 • • 



and, 



hi 


656 




■ 651 " 


_ 674 


676 




b^l 



'34 "36 



I2 

Scliur($i,$2) 
This concludes the proof of Corollarv l7.ll 



7.3. Singular examples. We now show how Theorem 15 . II allows to solve singular 
cases in Lidskii's theorem (Corollary 17. and we also illustrate the limitations of 
Theorem 01 

Example 7.3. Consider the following classical degenerate example, taken 
from |Wil65l Section 2.22] and |MB097I Eqn 1.1]: 



(61) Ae^Ao + eh, where Ao = N[];1] = 



■ 1 








1 




• 




□ • 






• 1 


□ • 




□ • 



and b G 



Recall that all the dots (whether they are surrounded by boxes or circles, or not) 
represent 0. If the entry &31 corresponding to the circled position in l|61(l . is zero, 
$1 is singular, and we cannot apply Lidskii's theorem (CoroUarv I7.1|l . However, 
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Theorem 15 . II can be applied. We can write {Ae)i 



with 





' bn 


1 


&13 


6i4 


615 " 




1 





1 


1 


1 




b2i 


b22 


1 


&24 


^25 




1 


1 





1 


1 


a = 





bz2 


^33 


^34 


635 


, and A = 


00 


1 


1 


1 


1 




bii 


bi2 


^43 


644 


1 




1 


1 


1 


1 







&51 


&52 


&53 


&54 


655 




1 


1 


1 


1 


1 



We have ai ~ Pm.in{A) = 2/5, Ci = {1,2,3,4,5}, and since the critical graph of 
A, which is composed only of the circuit (1^2^3^4^5-^1) covers all the 
nodes, we have G"=(i) = G"=(A). Thus, for G = G^(i), 





1 • 






• 1 








&45 ■ 






1 




b^i ■ ■ 





Theorem 15 . II shows that, if 645651 ^ Q, has five roots with asymptotics 

~ Ce'/' , where - 645651 • 



The asymptotics of the eigenvectors can also be obtained from Theorem 16.11 (the 
computations are similar to the case of Example 16. 2|) . 

Example 7.4. Let us discuss the following singular version of the illustrating ex- 
ample of |MB097| . Let At = ^0 + e6, where ^0 = -^[s^'i]' that, setting 
G = G'=(i), 



1 • 
• 1 

631 


634 • 


637 ■ 


639 


651 


664 • 


667 • 


669 


681 


684 • 


687 • 


689 


. 691 • 

case where 65 


694 • 
1 = 664 = 0. ^ 


697 • 

A'e may k 


699 . 
eep A 



but this gives little information since is not invertible. However. {Af)ij ~ aije^^^ 
still holds if we change the following values of A: A(,i = ^64 = 00. Then, we find 
ai = 1/3, Gi = {1,2,3}, a2 = 2/5, G2 = {4,5,6,7,8}, 03 = 4/5, G3 = {9}, and 
the critical graphs G^(A), t = 1,2,3 are represented as follows: 
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with the same colouring convention as in Example 14.51 

The matrix is invertible if, and only if, 631 7^ 0. In this case, Ae has three 
eigenvalues with asymptotics jC^ ~ Xe^^^, corresponding to the different cubic roots 
A of 631. We have 



" • 1 ■ 






■ 1 








1 










.Ki ■ ■ 




^99 . 



where for instance fog^ = 634 ^ ^8i^3i^&34- Thus, is invertible, if, and only if, 
^84^67 7^ 0. When this is the case, Ae has five eigenvalues with asymptotics ~ 
Ae^/^, corresponding to the different quintic roots A of 684667. 

The last critical graph, G'^{A), that we just represented above, does not have a 
disjoint circuit cover. To see this, observe that there is no arc from the set {3, 8, 9} 
to the set {2,3,5,6,7,8,9}, remark that the sum of the numbers of elements of 
these two sets, which is 3 + 7 = 10, exceeds the dimension of the matrix, which 
is 9, and apply the Frobenius-Konig theorem (see for instance |BR97I Th. 2.14]). 
Then, we know from Theorems l4.6l and l4.7l that the greatest root, 79, of the min-plus 
characteristic polynomial of A, Pa, is strictly greater than the greatest critical value, 
/3g = as = 4/5, and by Theorem 13.81 the exponent Ag of the remaining eigenvalue 
of Ae must be strictly greater than 4/5. Thus, in this case. Theorem 15 . II docs not 
predict the exponent of the eigenvalue of Ae of minimal modulus. However, this 
exponent can be obtained as follows. We already know that ai = 1/3 and 012 = 2/5 
are roots of respective multiplicity 3 and 5 of Pa, so the associated characteristic 
polynomial function is of the form pAiu) = [y ® cti)^{y ® "2)^(2/ © 79)- One can 
check that per^ = 4, and since Pa((D) = per A, we deduce that a\a\^'^ ~ 379 — 4, 
therefore, 79 = 1. Then, one can derive from Theorem 13.81 that A9 = 1, for generic 
values of b. The problem of finding the leading coefficient of the corresponding 
eigenvalue of Ae is solved by the result of |ABG04j . 

Example 7.5. Corollary 3.3 of |MB097j identifies a special situation where the lead- 
ing exponent of a group of eigenvalues can be found although the corresponding 
matrix <i>£_i appearing in Lidskii's theorem (see Corollarv l7. II above ^ is not invert- 
ible. We next give an example which cannot be solved using the method of |MB097j 
but which is solved by Theorem 15. II Let 



1 


6623 








£631 • 




£634 ■ 










1 


e656 












£^67 


£671 • 




£674 • 







where the bij are complex numbers, and, as above, the dots represent zero entries. 
The matrix Ae is of the form ^0 + d>, where the matrix Aq is a nilpotent matrix 
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conjugate to -^[21] ■ The corresponding matrix A is 



00 





00 


00 


00 


00 


00 


00 


00 
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00 


00 


00 


00 
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00 


00 
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00 


00 


00 


00 


00 


c» 


00 





00 


00 


00 


00 


00 


00 


00 


1 


00 


CO 


00 


00 


00 


00 


00 


1 


1 


00 


00 


1 


00 


00 


00 



We have ai = 2/3, Ci = {1,2,3}, 02 = 3/4, C2 = {4,5,6,7}. The critical graphs 
G1{A) and G2(A) are the following: 

1 1 




with the same colouring convention as in the previous examples. Theorem l5.1l shows 
that, for generic values of the coefficients bij, the matrix At has three eigenvalues 
£e ^ Ae^/'^, where A is a cubic root of 623631, and four eigenvalues of the form Ae^/^, 
where A is a quartic root of 

The following picture represents the actual Newton polygon of the characteristic 
polynomial det{YI — At), for generic values of the bij (this is exactly the graph of P, 
where P = per(Y/©y4.).) A monomial Y*e-' is represented by the point of coordinates 
(i, j). Integer points are represented by small crosses. The actual Newton polygon 
(black broken line) consists of two segments of respective slopes —3/4 and —2/3, 
joining the three circles. The approximation of the Newton polygon provided by 
Lidskii's theorem is given by the dashed broken line. 




The method of |MBU97j relies on the observation that Lidskii's theorem provides 
an approximation of the Newton polygon, which is exact when the matrices $^ 
arc mver tiblc. Corollary 3.3 of jMB097j requires the absence of integers points 
strictly between Lidskii's approximation and its chord, i.e. in the present case, in 
the interior of the gray region. Since this interior contains the integer point (4, 2), 
the leading exponents 2/3 and 3/4 cannot be obtained from |MB097) . 
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